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Abstract:  This  report  presents  results  of  a  three-dimensional  (3-D)  finite 
element  analysis  of  a  generic  miter  gate  with  multiple  cracks  used  to 
calculate  the  stress  intensity  factor  using  the  J-integral  approach.  The 
model  was  also  used  to  determine  the  optimum  diagonals  prestressing. 
Stress  intensity  factors  and  J-integral  approach  were  validated  numer¬ 
ically  against  closed-form  solutions  of  semi-infinite  plates  with  through¬ 
thickness  cracks,  double-edged  cracks,  single-edged  crack,  and  cracks 
growing  from  round  holes.  Only  a  2  percent  difference  between  the  closed- 
form  solution  and  the  numerical  calculated  values  was  obtained.  Addi¬ 
tionally,  the  report  shows  meshing  techniques  that  facilitated  the 
construction  of  the  complex  3-D  model.  Finally,  results  and  graphical 
methods  of  the  parametric  analysis  to  calculate  the  optimum  diagonal 
prestresses  are  illustrated. 
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1  Introduction 

1.1  Introduction 

Steel  hydraulic  structures  maintained  by  the  U.S.  Army  Corps  of  Engineers 
consist  primarily  of  lock  and  spillway  gates,  bulkheads,  and  other  closure 
structures.  In  recent  history  there  has  been  ample  evidence  of  unsatis¬ 
factory  performance  of  these  structures.  This  evidence  has  included 
numerous  cases  of  fatigue  cracking  caused  by  welded  connections  with  low 
fatigue  resistance  (Markland  Lock,  Greenup  Lock  and  Dam),  poor  weld 
quality,  unanticipated  structural  behavior,  or  unexpected  loading.  In  the 
most  severe  cases  at  least  one  catastrophic  failure  of  a  bulkhead  system 
nearly  resulted  in  loss  of  life  (Coffeeville  Lock,  maintenance  stoplogs, 
1994),  and  complete  replacement  of  lock  gates  has  been  required  in  at 
least  two  other  cases  (Melvin  Price  Lock  and  Upper  Mississippi  River 
Lock  27).  Miter  gates  have  also  been  investigated  and  cracks  repaired. 
Maintenance  and  repair  of  fatigue  and  fracture  failures  represents  a  major 
Operations  and  Maintenance  expenditure  for  the  Corps. 

Because  of  these  failures,  analytical  techniques  for  employing  state-of-the- 
art  capabilities  for  fracture  mechanics  analysis  using  finite  element 
modeling  have  been  developed.  This  effort  included  the  use  of  commer¬ 
cially  available  nonlinear  finite  element  programs  incorporating  J-integral 
analysis  for  fracture  analysis.  The  results  of  J-integral  analysis  can  be 
directly  compared  to  elastic  or  elastic-plastic  material  properties  for 
assessment.  These  analytical  tools  can  be  used  for  much  more  accurate 
and  detailed  fracture  assessments  than  the  Corps  currently  has  the 
capability  to  perform  and  can  be  used  to  access  the  criticality  or  need  to 
perform  repairs. 

Detailed  finite  element  fracture  analyses  employing  J-integral  analysis 
have  the  potential  to  provide  a  much  more  accurate  fracture  assessment. 
Cost  savings  may  be  realized  from  accurate  fracture  assessments  through 
avoidance  of  unnecessary  repairs.  Detailed  fracture  analysis  can  also 
provide  an  important  investigative  tool  in  assessing  the  actual  loads  in  a 
member  at  the  time  of  failure. 
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Recently,  there  have  also  been  a  number  of  failures  of  diagonals  on  miter 
gates  (Figure  l).  These  appear  to  be  fatigue-induced  failures  driven  by  the 
connection  details  at  the  diagonal  ends,  but  some  evidence  also  suggests 
that  the  current  design  guidance  results  in  a  much  larger  prestress  in  miter 
gate  diagonals  than  may  be  required.  The  high  prestress  theoretically 
increases  the  fatigue  of  these  members  and  results  in  increased  stress  in 
the  miter  gate. 


Figure  1.  Diagonal  failure. 


Comprehensive  three-dimensional  (3-D)  finite  element  analyses  have  the 
capability  to  assess  the  diagonal  prestressing  levels,  which  will  provide  the 
initiation  for  a  new  design  approach.  The  new  techniques  will  allow  the 
prediction  of  optimum  prestressing  levels  without  compromising  the 
fatigue  capacity. 

Most  recently,  a  number  of  gates  have  experienced  cracking  problems 
around  the  pintle  region  (Figure  2),  in  addition  to  buckling  occurring  in 
vertical  plates  on  the  miter  end  sections.  The  3-D  finite  element  modeling 
will  provide  the  basis  to  understand  the  source  of  the  problem  and 
generate  appropriate  repair  procedures. 

1.2  Objectives 

The  proposed  effort  includes  the  development  of  analytical  techniques  for 
employing  state-of-the-art  capabilities  for  fracture  mechanics  analysis 
using  finite  element  modeling.  This  effort  includes  the  use  of  commercially 
available  nonlinear  finite  element  programs  incorporating  J-integral 
analysis  for  fracture  analysis. 


ERDC/ITL  TR-09-1 


3 


Figure  2.  Cracking  in  pintle  zone. 


The  results  of  J-integral  analysis  can  be  directly  compared  to  elastic  or 
elastic-plastic  material  properties  for  assessment.  These  analytical  tools 
can  be  used  for  much  more  accurate  and  detailed  fracture  assessments 
than  the  Corps  currently  has  the  capability  to  perform.  An  effort  to  assess 
the  connection  detailing  and  design  of  miter  gate  diagonals  is  proposed. 
This  effort  will  include  development  of  analytical  models  to  assess  the 
prestress  requirements  of  diagonals  and  to  develop  improved  details  for 
diagonal/gate  connections.  Last,  the  pintle  region  cracking  and  buckling 
problem  will  be  examined,  and  recommendations  for  predicting  the 
problem  will  be  discussed. 

1.3  Research  significance 

The  Corps  requires  the  inspection  of  hydraulic  steel  structures  with  frac¬ 
ture  critical  members  (FCM)  whose  failure  would  result  in  probable  loss  of 
life  to  ensure  structural  integrity.  Detailed  finite  element  fracture  analyses 
employing  J-integral  analysis  have  the  potential  to  provide  a  much  more 
accurate  fracture  assessment.  The  benefits  of  an  accurate  fracture  assess¬ 
ment  may  be  realized  from  cost  savings  through  avoidance  of  unnecessary 
repairs.  Detailed  fracture  analysis  can  also  provide  an  important  investi¬ 
gative  tool  in  assessing  the  actual  loads  in  a  member  at  the  time  of  failure. 
The  reassessment  of  miter  gate  diagonal  design  criteria  and  connection 
details  should  lead  to  more  realistic  prestress  levels  in  diagonals  and  could 
provide  the  initiation  of  new  design  criteria  that  eliminate  the  potential  of 
fatigue  failures  in  the  diagonal  connection.  Considering  the  number  of 
miter  gates  exhibiting  pintle  region  cracking  and  buckling,  understanding 
the  causes  of  the  problem  is  not  just  necessary  but  essential. 
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1.4  Overview 

Chapter  2  includes  a  brief  review  of  the  basic  concepts  of  linear  elastic 
fracture  mechanics.  It  also  presents  a  detailed  explanation  of  the  J-integral 
approach.  Chapter  3  presents  the  validation  of  the  J-integral  computed 
numerically  against  closed-form  solutions  of  semi-infinite  plates  with 
through-thickness  cracks,  double-edged  cracks,  single-edged  cracks,  and 
cracks  growing  from  round  holes.  Details  of  the  numerical  models  devel¬ 
oped  in  this  investigation  are  presented  in  Chapter  4.  Results  from  this 
study  along  with  analysis  of  these  results  are  presented  in  Chapter  5. 
Conclusions  and  recommendations  for  further  studies  are  presented  in 
Chapter  6. 
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2  Fracture  Mechanics  and  its  Applications 
to  Finite  Elements  in  Steel  Structures 

2.1  Introduction 

Fracture  mechanics  is  a  method  of  characterizing  the  fracture  behavior  in 
structural  parameters,  stress,  and  flaw  size,  which  can  be  used  directly  as 
Barsom  and  Rolfe  stated  (1987).  They  also  declared  that  the  science  of 
fracture  mechanics  can  be  used  to  describe  quantitatively  the  trade-offs 
among  stresses,  material  toughness,  and  flaw  size. 

The  concepts  of  fracture  mechanics  were  used  in  the  World  War  II  era  to 
predict  the  failure  of  brittle  steel  materials  in  boat  hulls.  Since  then,  linear 
elastic  fracture  mechanics  (LEFM)  has  been  increasingly  used  to  predict 
brittle  fracture  of  high-strength  low-toughness  materials. 

The  main  purpose  is  to  determine  the  rates  of  change  of  an  existing  crack, 
if  this  crack  will  grow  at  a  certain  load  magnitude,  and  if  it  does,  at  what 
rate  and  direction.  The  numerical  simulation  requires  the  state  of  stresses, 
strain,  and  displacements  present  around  the  crack  tip  from  which  it  can 
be  predicted  if  the  crack  will  grow  or  not,  and  if  it  grows,  in  what  direction. 
For  LEFM  and  the  work  described  herein,  the  numerical  simulation  will  be 
used  to  obtain  the  stress  intensity  factor  for  Mode  I.  This  chapter  will 
describe  the  basic  concepts  of  LEFM  and  the  method  for  obtaining  the 
stress  intensity  factor  for  Mode  I  with  the  J-integral  approach. 

2.1  Modes  of  Fracture 

There  are  three  principal  modes  of  fracture,  commonly  referred  to  as  the 
opening  mode,  the  sliding  mode,  and  the  tearing  mode.  The  difference 
between  these  modes  of  fracture  is  in  the  type  of  relative  displacement 
between  the  surfaces  of  a  crack. 

2.2.1  Mode  I,  the  opening  mode 

In  the  opening  mode,  the  instantaneous  displacement  of  two  crack  sur¬ 
faces  is  in  opposite  directions  and  perpendicular  to  the  crack  surface  as 
shown  in  Figure  3a.  Direct  tension  and  flexural  action  usually  cause 
structural  failure  in  this  mode. 
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Figure  3.  Fracture  failure  modes,  (a)  Mode  I  -  opening  mode;  (b)  Mode  II  -  sliding 
or  shear  mode;  (c)  Mode  III  -  tearing  or  out-of-plane  shearing  mode. 

2.2.2  Mode  II,  the  sliding  or  shearing  mode 

In  the  sliding  or  shearing  mode,  instantaneous  movement  of  the  crack 
surfaces  is  in  the  opposite  direction  and  perpendicular  to  the  crack  front 
line  but  in  the  plane  of  the  crack  surface.  So,  plane  x-y  is  plane  of  sym¬ 
metry,  but  cracking  with  respect  to  plane  x-z  is  antisymmetric.  This  mode 
of  fracture  is  usually  observed  when  in-plane  shear  force  is  imposed  on  the 
structure  (Figure  3b). 

2.2.3  Mode  III,  the  tearing  mode 

Figure  3c  shows  the  tearing  mode.  In  this  mode,  instantaneous  move¬ 
ments  of  the  two  crack  surfaces  are  in  opposite  directions  but  in  the 
direction  of  the  crack  front  line,  and  in  the  plane  of  crack  surface,  it  is 
antisym-metric  with  respect  to  both  planes  x-y  and  x-z.  Torsion  loading 
causes  this  mode  of  fracture. 

2.3  Linear  Elastic  Fracture  Mechanics  (LEFM) 

Linear-elastic  fracture-mechanics  technology  is  based  on  an  analytical 
procedure  that  relates  the  stress-field  magnitude  and  distribution  in  the 
vicinity  of  a  crack  tip  to  the  nominal  stress  applied  to  the  structural 
member;  to  the  size,  shape,  and  orientation  of  the  crack  or  crack-like 
discontinuity;  and  to  the  material  properties  (Barsom  and  Rolf  1987). 
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The  fundamental  principle  of  fracture  mechanics  is  that  the  stress  field 
ahead  of  a  sharp  crack  in  a  structural  member  can  be  characterized  in 
terms  of  a  single  parameter,  K,  the  stress  intensity  factor  that  has  units 
ofpsiViir  This  parameter,  K,  is  related  to  both  the  nominal  stress  level  o  in 
the  member  and  the  size  of  the  crack  present  a. 


LEFM  is  the  basic  theory  of  fracture,  originated  by  Griffith  (1921)  and 
completed  by  Irwin  (1957)  and  Rice  (1968a,  1968b).  The  LEFM  is  a  highly 
simplified  theory  that  is  applicable  to  any  material  as  long  as  certain  con¬ 
ditions  are  met.  These  conditions  are  related  to  the  basic  ideal  situation 
analyzed  in  LEFM  in  which  all  the  material  is  elastic  except  in  a  vanish¬ 
ingly  small  region  (a  point)  at  the  crack  tip  (Bazant  and  Planas  1998). 


2.3.1  Griffith’s  fracture  criterion 

If  the  amount  of  energy  supplied  to  the  crack  tip  for  an  incremental  crack 
growth  is  equal  to  or  greater  than  the  amount  of  energy  required  for  crack 
growth,  then  this  results  in  catastrophic  crack  growth  (Anderson  1991). 

The  concept  forms  the  basis  of  LEFM.  The  energy  supplied  to  the  crack  tip, 
U,  is  usually  either  the  energy  released  from  crack  formation  or  the  work 
done  by  the  applied  forces  from  additional  displacement  caused  by 
introducing  a  crack.  W  denotes  the  amount  of  energy  absorbed  in  the 
creation  of  new  free  surfaces  as  a  result  of  crack  propagation.  The  Griffith 
criterion  for  catastrophic  crack  growth  is  then  given  by: 


dU_>dW_ 
da  da 


(2.1) 


where  da  is  the  incremental  crack  growth  and  dU /da  is  the  energy  release 
rate,  G.  For  a  center  crack  plate  (crack  length  2 a)  subjected  to  far-field 
tensile  stress,  o,  it  can  be  shown  that: 


dU  _  a2jia 
da  E 


(2.2) 


where  E  is  the  modulus  of  elasticity,  and  dW/da  is  crack  resistance,  R, 
given  by: 


R 


dW 


2  y 


da 


(2.3) 
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where  y  is  the  specific  surface  energy,  considered  a  material  property. 
Consequently,  if 


a2Jia 

E 


>2  y 


or 


cr> 


2Ey 
ji a 


the  crack  will  grow  catastrophically. 

The  energy  release  rate,  G,  for  Mode  I  fracture  is  denoted  as  Gi.  At 
incipient  catastrophic  crack  growth,  Gi  equals  Gic,  which  is  the  critical 
energy  release  rate  (Gic  =  2y).  Gic  is  considered  a  material  property. 

2.3.2  Stress  distribution  ahead  of  the  crack  tip 

In  1898,  Kirsch  obtained  the  stress  concentration  at  the  edge  of  a  circular 
hole  as  three  times  the  nominal  stress  (Figure  4a).  In  1913,  the  same  prob¬ 
lem  was  solved  for  a  more  general  case  of  an  elliptical  hole  by  Ingliss 
(Inglis  1913)  (Figure  4b).  For  that  case,  the  stress  concentration  factor  was 
obtained  as  1  +  2  *  a/b .  A  sharp  crack  in  a  solid  body  can  be  approximated 

in  the  limiting  case  with  an  elliptical  crack  where  b  — >  o  (Figure  4c).  The 
related  expression  for  the  stress  distribution  shows  that  the  normal  stress 
at  the  crack  tip  will  be  infinity  for  a  line  crack. 
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Figure  4.  Stress  concentration  at  the  crack  tip  due  to  (a)  circular  hole,  (b)  elliptical  hole,  and 

(c)  sharp  line  crack. 


Irwin  (Irwin  1957)  was  credited  with  obtaining  mathematical  expressions 
for  the  stress  field  near  the  crack  tip.  Stress  field  (a,))  at  the  crack  tip  in 
each  of  the  principal  modes  of  fracture  can  be  obtained  in  the  form: 
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K 


jir 


(2.4) 


where  r  is  the  distance  from  the  crack  tip  and  cp  is  the  angle  of  the  stress 
tensor yij(cp),  as  shown  in  Figure  5.  Equation  2.4  is  valid  for  small  values  of 
r.  The  form  of fij  depends  on  the  specimen  and  crack  geometry.  K  in 
Equation  2.4  is  called  the  stress  intensity  factor,  and  equals  to  Ki,  Kn,  or 
Km  depending  on  the  applicable  mode  of  fracture  (I,  II,  or  III, 
respectively). 


t  t  t  t  t  t  t  t  T  t  t  t  t  t  t  t  t 


Figure  5.  Notations  for  stresses  and  other  geometric  parameters. 

Irwin  found  that  the  stress  and  displacement  fields  in  the  vicinity  of  crack 
tips  subjected  to  the  three  modes  of  deformation  are  given  by  Equations 
2.5,  2.6,  and  2.7.  Irwin  developed  these  equations  using  the  Westergaard 
method.  Figure  4  shows  the  stress  concentration  at  the  crack  tip  caused  by 
(a)  a  circular  hole,  (b)  an  elliptical  hole,  and  (c)  a  sharp  line  crack.  The 
stress  distribution  for  Mode  I  fracture  (Barsom  and  Rolfe  1987),  can  be 
obtained  as: 


V^rr 


V 

1  _ qiri 

cp) 

q-iri 

(3cp j 

2 , 

_L  Olll 

Uj 

Olll 

2  J] 

(2.5) 
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G 


Kt 


yy 


>/2 


■  cos 


jrr 


<P 

\  2  / 


1  +  sin 


( q)^ 


sin 


{3cp ^ 


(2.6) 


Txy  V2 


r  COS 


nr 


sin 


sm 


(2.7) 


At  cp  =  o  one  obtains 


cr 


Kr 


yy 


*?i*r 


The  distribution  obtained  for  ayy  is  also  shown  in  Figure  5. 

2.3.3  Stress  intensity  factor 

The  second  major  achievement  in  the  theoretical  foundation  of  LEFM  was 
the  introduction  of  the  stress  intensity  factor,  K,  as  a  parameter  for  the 
intensity  of  stresses  close  to  the  crack  tip  and  related  to  the  energy  release 
rate  (Bazant  and  Planas  1998).  Ingliss  (1913)  studied  the  unexpected  fail¬ 
ure  of  naval  ships,  and  Griffith  (1921)  extended  this  work  using  thermo¬ 
dynamical  criteria.  Using  this  work  Irwin  (1957)  developed  this  concept  of 
the  stress  intensity  factor. 

The  stress  intensity  factor  is  a  measure  of  the  change  in  stress  within  the 
vicinity  of  the  crack  tip.  Therefore,  it  is  important  to  know  the  crack 
direction  and  when  the  crack  stops  propagating.  The  stress  intensity  factor 
is  compared  with  the  critical  stress  intensity  factor  to  determine  whether 
or  not  the  crack  will  propagate. 

Dimensional  analysis  can  be  used  to  show  that  the  stress  intensity  factor 
for  Mode  I  fracture  Ki  has  the  following  form: 


Kj  —  gO'Jjm  (2.8) 

where  g  is  a  nondimensional  function  depending  on  the  size  and  geometry 
of  the  crack,  size  and  geometry  of  the  structural  component,  and  the  type 
of  loading.  For  normal  cracks,  its  value  ranges  between  1  and  2,  but  maybe 
larger  for  longer  cracks.  Functions  defined  for  common  geometries  and 
loading  conditions  are  available  in  published  literature  (Barsom  and  Rolfe 
1987;  Tada  et  al.  1985).  a  is  the  nominal  far-field  stress  and  2a  is  the  crack 
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length.  If  Ki  is  the  same  for  two  cracked  bodies,  then  based  on  the  equa¬ 
tions,  the  same  stress  field  will  exist  at  their  crack  tips.  If  the  two  bodies 
are  made  of  the  same  material,  an  identical  response  is  expected.  This  fact 
results  in  the  important  conclusion  that  Ki  can  be  used  as  a  similitude 
parameter  to  compare  the  response  of  the  same  material  at  the  crack  tip 
and  also  to  compare  the  degree  to  which  materials  are  influenced  by  the 
stress  fields. 

2.3.4  Fracture  toughness:  Critical  stress  intensity  factor 

Another  important  parameter  of  the  LEFM  is  the  fracture  toughness,  Kic 
(the  capacity).  Fracture  toughness  is  a  material  fracture  property  (Swamy 
1979).  If  a  plate  is  loaded  (Figure  6)  to  the  failure  stress,  a/,  and  if  a  =  <7/ is 
the  value  of  the  nominal  far-field  stress  at  failure,  then  a  value  of  Ki  associ¬ 
ated  with  Of  could  be  determined  and  referred  to  as  Kic.  If  this  toughness 
property  is  available  for  the  material,  then  the  failure  state  stress  could  be 
expressed  as  (Swamy  1979): 


a 


f  ~ 


(2.9) 


ttttttt 


TTTTTT 


Figure  6.  Failure  stress  related  to  crack  size  for  infinitely  wide  plate  subjected  to  tension. 


From  Figure  6  it  can  be  observed  that  when  the  obtained  value  of  07 
exceeds  the  tensile  strength  of  the  material,  ft  (or  when  a  <  ac,  the  critical 
crack  length),  for  small  crack  length,  the  LEFM  is  not  valid.  The  relation 
between  the  energy  release  rate,  Gi,  and  the  stress  intensity  factor,  Ki,  can 
be  shown  (Barsom  and  Rolfe  1987;  Bocca  et  al.  1991)  to  be  of  the  form: 
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Gj=^h  (2.10) 

hi 

where  for  plane  stress  conditions  E'  =  E  and  for  plane  strain  conditions 
E'=E/(i-v2).  Where  E  is  the  modulus  of  elasticity  and  v  is  the  Poisson 
ratio.  Instead  of  writing  the  catastrophic  crack  growth  condition  in  terms 
of  Gi  (as  Gi  =  Gic )  one  may  alternatively  use  the  condition  Ki  =  Kic  ( Kic  can 
be  obtained  using  Equation  2.10  by  substituting  Gi  =  G/c).  Kic,  and  G/c  are 
not  loading  and  crack  geometry  dependent.  Like  G/c,  Kic,  also  termed  the 
fracture  toughness  of  the  material,  is  a  material  property. 

2.3.5  Plastic  zone  and  influence  of  constraint  at  the  crack  tip 

Realistically,  using  elastic  theories,  the  materials  cannot  be  expected  to 
sustain  the  extremely  large  stresses  predicted  at  the  crack  tip.  The  theo¬ 
retical  shape  of  the  local  yield  (plastic)  zone  can  be  determined  using  the 
elastic  stress  distribution  obtained  in  the  vicinity  of  the  crack  tip  and  the 
yield  strength  of  metals  (shaded  zone  in  Figure  5).  Elongation  in  the 
direction  of  applied  stress  is  accompanied  by  negative  strains  in  the  other 
perpendicular  directions  from  the  influence  of  the  Poisson’s  ratio.  In  other 
words,  an  elongation  occurs  in  the  applied  stress  direction  and  shrinkage 
in  the  other  direction.  Plane  stress  or  plane  strain  conditions  may  exist  at 
the  crack  tip,  depending  on  the  type  of  constraint  at  the  crack  tip. 

Usually  at  the  surface,  the  plane  stress  condition  results  in  a  large  plastic 
zone  from  the  low  level  of  constraint.  Also,  with  the  large  plastic  defor¬ 
mation  and  more  energy  absorption,  there  is  more  resistance  to  crack 
growth.  Plane  strain  conditions  prevail  in  the  interior  of  the  specimen, 
where  constraints  are  more  severe.  A  smaller  plastic  zone  with  less  energy 
absorption  is  typical  under  plane  strain  conditions  of  crack  growth.  The 
thickness  of  the  specimen,  as  a  result,  plays  a  major  role  in  determining 
the  state  of  stress  at  the  crack  tip. 

2.4  J-integral 

The  J-integral  was  presented  by  Cherepanov  (1967)  and  by  Rice  (1968a, 
1968b)  for  two-dimensional  (2-D)  domains  containing  cracks.  The 
J-integral  is  usually  used  in  rate-independent  quasi-static  fracture  analysis 
to  characterize  the  energy  release  associated  with  crack  advance.  It  can  be 
related  to  the  stress  intensity  factor  if  the  material  response  is  linear. 


ERDC/ITL  TR-09-1 


14 


Consider  a  2-D  linear  body  of  linear  or  nonlinear  elastic  material  free  of 
body  forces  and  subjected  to  a  2-D  deformation  field  (plane  strain,  plane 
stress)  so  that  all  stresses  Oij  depend  only  on  two  Cartesian  coordinates 
(x,  y).  Suppose  the  body  contains  an  edge  crack  as  shown  in  Figure  7. 


Figure  7.  Crack  tip  coordinate  system  and  typical  line 
integral  contour. 


Define  the  strain-energy  density  W  as 


(2.11) 


where  e  =  [e,j]  is  the  infinitesimal  strain  tensor.  Now,  define  the  J-integral 
as 


(2.12) 


As  shown  in  Figure  7,  r  is  the  curve  surrounding  the  notch  tip,  the  integral 
being  evaluated  counterclockwise  starting  from  the  lower  surface  and 
continuing  along  the  path  r  to  the  upper  surface.  In  Equation  2.12,  T  is  the 
traction  vector  defined  according  to  the  outward  normal  along  r, 

X)  =  o^Tij :  u  is  the  displacement  vector,  and  ds  is  an  element  arc  length 

along  r.  Rice  (1968a,  1968b)  proved  the  path-independent  concepts  and 
found  that  for  small-scale  yielding  the  stress  energy  release  rate  G  is  equal 
to  the  J-integral.  Therefore,  the  stress  intensity  factor  can  be  evaluated  as 
follows: 
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G  =  ^r 

E 

(2.13) 

K2, 

J-G-  { 

E 

(2.14) 

[| 

II 

(2.15) 

where: 

E  =  -  —  -  for  plane  strain  and  E*  =  E  for  plane  stress 


For  a  virtual  crack  advance  cp(p)  in  the  plane  of  a  three-dimensional 
fracture,  the  energy  release  rate  is  given  by 

J  =  f  ( p)n-Q  rdA  (2.16) 

A 


where  dA  is  a  surface  element  along  a  vanishing  small  tubular  surface 
enclosing  the  crack  tip  or  crack  line,  n  is  the  outward  normal  to  dA,  and  r 
is  the  local  direction  of  virtual  crack  extension.  Q  is  given  by 


Q  = 


1  du  ^ 

WI-o  — 
dx 


(2.17) 


For  elastic  material  behavior  W  is  the  elastic  strain  energy;  for  elastic- 
plastic  or  elasto-viscoplastic  material  behavior  W  is  defined  as  the  elastic 
strain  energy  density  plus  the  plastic  dissipation,  thus  representing  the 
strain  energy  in  an  equivalent  elastic  material.  Therefore,  the  J-integral 
calculated  is  suitable  only  for  monotonic  loading  of  elastic-plastic 
materials.  I  is  the  unit  matrix. 


The  stress  intensity  factors  Ki,  Kn,  and  Km  are  usually  used  in  linear  elastic 
fracture  mechanics  to  characterize  the  local  crack-tip/ crack-line  stress  and 
displacement  fields.  They  are  related  to  the  energy  release  rate  (the 
J-integral)  through 


J  =  —KtB~1K 
8jt 


(2-18) 
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where  K  =  [ Ki ,  Ku,  Km]  are  the  stress  intensity  factors  and  B  is  called  the 
prelogarithmic  energy  factor  matrix.  For  homogeneous,  isotropic  materials 
B  is  diagonal,  and  Equation  2-18  simplifies  to 

J  =  ±(Kl+Kl)+±Kl„  (2-19) 

where  E  =  E  for  plane  stress  and  E  =  E/(  1-  v 2)  for  plane  strain, 
axisymmetry,  and  three  dimensions. 
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3  Fracture  Mechanics  Numerical  Solution: 
Validation  of  J-integral  Approach 

The  field  of  fracture  mechanics  has  become  the  primary  approach  to 
control  brittle  fracture  and  fatigue  failures  in  structures.  An  overwhelming 
amount  of  research  on  brittle  fracture  in  structures  of  all  types  has  shown 
that  numerous  factors  can  contribute  to  brittle  fractures  in  large  welded 
structures,  for  example,  service  temperature,  material  toughness,  connec¬ 
tions  details,  welding  residual  stresses,  fatigue,  and  constraint.  However 
some  primary  factors  control  the  susceptibility  of  a  structure  to  brittle 
fracture:  the  material  toughness,  Kic,  the  crack  size  ac;  and  the  stress 
level  a. 

The  material  toughness  ( Kic )  can  be  defined  as  the  ability  to  carry  load  or 
deform  plastically  in  the  presence  of  a  notch  and  can  be  described  in  terms 
of  the  critical  stress  intensity  factor  under  conditions  of  plane  stress,  plain 
strain,  or  critical  toughness. 

The  critical  stress  intensity  factor,  Kic,  represents  the  terminal  conditions 
in  the  life  of  a  structural  component.  The  total  useful  life  of  a  component  is 
determined  by  the  time  necessary  for  a  crack  to  initiate  and  propagate 
from  subcritical  dimensions  to  the  critical  size,  ac. 

This  chapter  presents  the  validation  of  the  J-integral  computed  numer¬ 
ically  against  closed-form  solutions  of  semi-infinite  plates  with  through¬ 
thickness  cracks,  double-edge  cracks,  single-edge  cracks,  and  cracks 
growing  from  round  holes. 

3.1  Through-thickness  crack 

A  perfectly  rectangular  plate  of  20  by  40  in.  (Figure  8)  with  a  center  crack 
of  5.0  in.  with  a  tensile  stress  of  100.0  psi  was  analyzed  and  compared 
with  the  closed  form  solution  (Tada  et  al.  1985).  Figure  9  shows  the  finite 
element  mesh  and  stress  contours  for  the  finite  element  analysis. 
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Figure  8.  Plate  sketch  and  dimensions  for  through-thickness  crack. 


Figure  9.  Through-thickness  crack  plate  finite  element  mesh  and 
stress  contour. 


The  stress  intensity  factor  for  an  infinite  plate  subjected  to  uniform  tensile 
stress,  a,  containing  a  through-thickness  crack  of  length  2 a  is: 


Kj  =  OyJjia 


(3.1) 


ERDC/ITL  TR-09-1 


19 


with  a  tangent  correction  of 


26  ^  Jia 
—  tan — 

jia  2b 

\  / 


1/2 


(3.2) 


It  becomes 


Kj  =  o^jjra 


2b  ^  Jia 
—  tan— - 
jia  2b 


1/2 


(3.3) 


To  obtain  the  correction  factor,  an  interpolation  is  performed  using  the 
values  in  Table  l: 


a  =  2.5" 
6  =  10" 


a 

6 


2.5 


10 


-0.25 


2b  ^  Jia 
—  tan — 
jia  2b 


xl/2 


1.03 


Table  1.  Correction  factors  for  a  finite-width  plate  containing  a  through-thickness  crack. 


a/b 

( \i/2 

26  ^  Jia 
—  tan — 

jia  2b 

\  / 

a/b 

(  \i/2 

26  ^  Jia 
—  tan— - 

jia  2b 

\  / 

0.074 

1.00 

0.410 

1.08 

0.207 

1.02 

0.466 

1.11 

0.275 

1.03 

0.535 

1.15 

0.337 

1.05 

0.592 

1.20 

Then, 

Kj  =  (100)  ^71-  2.5 -(1.03) 

Kj  =  288.657 psiy/in 

The  results  from  the  finite  element  evaluation  provide  a  J-integral  value  of 
2.82  x  10  3  psi-in.  Using  Equation  2.14  to  determine  the  stress  intensity 
factor  Ki  gives: 
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E 

J  =  2.82  x  10'3  psi  -  in 
E  =  30  xlO6  psi 


K  =  ylJ-E  =  V2.82x10“3-30x106  « 290.86  psiVin. 

The  results  show  a  difference  of  less  than  l  percent  between  the  numerical 
and  closed-form  solutions. 

3.2  Double-edge-notched  plate 

A  perfect  rectangular  plate  of  20  by  40  in.  (Figure  10)  with  double  edge 
cracks  of  2.5  in.  with  a  tensile  stress  of  100.0  psi  was  analyzed  and  com¬ 
pared  with  the  closed-form  solution.  Figure  11  shows  the  finite  element 
mesh  and  stress  contours  for  the  finite  element  analysis. 


Figure  10.  Plate  sketch  and  dimensions  for  double-edge-notched  plate. 
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Figure  11.  Double-edge-notched  plate  finite  element  mesh  and 
stress  contour. 


The  stress  intensity  factor  for  double-edge-notched  infinite  plate  that  is 
subjected  to  a  uniform  tensile  stress  is  given  by k,  =1.12 ajna ,  where  the 

constant  1.12  represents  a  free  surface  correction  factor  for  edge  notches 
that  are  perpendicular  to  the  applied  tensile  stress.  For  this  example, 

a  =2.5  in. 

Kj  =  1.12  a^JJTa 
Kj  =1.12  (100)  72.5  ^ji 
Kj  =  313.88  psi^Jin. 

The  results  from  the  finite  element  evaluation  provide  a  J-integral  value  of 
3.26  x  iO"3  psi-in.  for  the  left  edge  crack  and  3.18  x  io-3  psi-in.  for  the  right 
edge  crack.  For  this  example  a  similar  value  of  the  J-integral  would  be 
expected  for  both  cracks;  however,  because  of  the  lack  of  mesh  symmetry 
the  values  differ.  It  is  important  to  point  out  that  the  J-integral  is  directly 
proportional  to  the  strain  energy  density  of  those  finite  elements  inside  the 
contour  path  chosen;  therefore,  the  lack  of  mesh  symmetry  will  vary  the 
resultant  J-integral  as  shown  in  this  example.  Then  using  Equation  2.14  to 
determine  the  stress  intensity  factor  Ki  produces  less  than  1  percent 
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difference  is  observed  between  the  mean  value  of  the  stress  intensity  factor 
obtained  from  the  finite  element  analysis  and  the  closed-form  solution. 


E 

J  =  3.26x10  ~3 psi -in 
£  =  30x10®  psi 


K  =  VjTe  =  ^3.26  x  10-3  •  30  x  106  «  312.73  psiVim  for  the  left  crack 
and: 


J  =  3.18x10  3psi-in 
E  =  30  xlO6  psi 


K  4T~E  73.18  10  3  30  106  308.87  psi  Vac  for  the  right  crack. 

3.3  Single-edge-notched  plate 

For  this  example,  a  plate  with  the  same  40  by  20  in.  dimensions  as  that  of 
the  previous  examples  is  used  (Figure  12).  A  single-edge  crack  of  2.5  in. 
was  induced  at  the  left  edge  of  the  plate  and  the  plate  was  loaded  with  100 
psi  in  tension  as  shown  in  Figure  12.  Figure  13  shows  the  finite  element 
mesh  and  stress  contours  for  the  finite  element  analysis. 

The  stress-intensity-factor  equation  for  a  semi-finite  single-edge  notched 
specimen  is  given  by  Equation  3.5: 


Kj=  1.12a  Jjra-k 


(3.5) 


This  function  has  an  additional  correction  factor  to  account  for  bending 
stresses  caused  by  lack  of  symmetry  in  the  single-edge-notched  specimen. 
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a 


Figure  12.  Plate  sketch  and  dimensions  for  single-edge-notched  plate. 


stress  contour. 


To  obtain  the  correction  factor,  an  interpolation  is  performed  using  the 
values  in  Table  2: 


a  =  2.5  in. 
b  =  10  in. 


1.11 


a 


2.5 


0.25 


b  10 
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Table  2.  Correction  factors  for  a  single-edge-notched  plate. 


a/b 

k(a/b) 

a/b 

k(a/b) 

0.10 

1.03 

0.60 

1.50 

0.20 

1.07 

0.70 

1.69 

0.30 

1.15 

0.80 

1.91 

0.40 

1.22 

0.90 

2.20 

0.50 

1.35 

1.00 

2.55 

Then, 

Kj  =1.12  (100)  y/n- 2.5 -(1.11) 

Kj  =  348.406  psiVinT 

The  results  from  the  finite  element  evaluation  provide  a  J-integral  value  of 
3.92  x  10-3  psi-in.  Using  Equation  2.14  to  determine  the  stress  intensity 
factor  Ki,  gives 

j-* 

E 

J  =  3.92  xl0“3  psi-in 
E  =  30xl06psi 


K  =  VjTe  =  ^3.92x10“3-30x106  «  342.93  psiVim 

Once  again,  a  difference  of  less  than  1  percent  is  observed  between  the 
stress  intensity  factor  from  the  finite  element  analysis  and  the  closed-form 
solution. 

3.4  Cracks  growing  from  round  holes 

For  this  example,  a  rectangular  plate  40  in.  long  (h  =  40  in.)  and  20  in. 
wide  ( b  =  20  in.)  with  a  5.o-in.-radius  hole  and  two  1.0-in.  edge  cracks 
extending  toward  the  longitudinal  edges  of  the  plate  were  introduced 
(Figure  14).  The  plate  was  loaded  with  100.0  psi  in  tension.  Figure  15 
shows  the  finite  element  mesh  and  stress  contours  for  the  finite  element 
analysis. 
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CT 


CT 


Figure  14.  Plate  sketch  and  dimensions  for  cracks  growing 
from  round  holes. 


Figure  15.  Crack  before  analysis  for  cracks  growing  from  round  holes. 


The  stress  intensity  factor  for  cracks  growing  from  a  circular  hole  in  an 
infinite  plate  is  given  by: 


(3.6) 
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Where  for  our  example, 


a  =  6  in. 

6  =  10  in. 

a  6  .  _ 

—  =  —  =  0.6 

6  10 


h  =  20  in. 

6  =  10  in. 

h  20  ^ 

-  =  —  =  2.0 

6  10 


R  =  5  in. 

6  =  10  in. 

a  5  _  _ 

—  =  —  =  0.5 

6  10 


a  = 100  psi 


The  value  for  F 


a  R  h 
b’~b’b 


is  obtained  from  Tada  et  al.  (1985)  as  1.5,  then, 


K  —  oJna  F 


a  R  h 
b’b’b 


1007^6  (1.5)  =  308.28  psiVIn. 


The  results  from  the  finite  element  evaluation  provide  a  J-integral  value  of 
3.05  x  10  3  psi-in.  Using  Equation  2.14  to  determine  the  stress  intensity 
factor  Ki,  gives 


E 


J  =  3.05  x  10  3  psi-in. 
E  =  30xl06psi 


K  =  VJ-E  =  73.05  xl0“3- 30 xlO6  ^  302.49  psiTin. 

For  this  example,  a  difference  of  less  than  2  percent  was  observed  between 
the  stress  intensity  factor  from  the  finite  element  analysis  and  the  closed- 
form  solution. 
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4  3-D  Finite  Element  Analysis  of  Miter  Gate 

With  Multiple  Cracks 

This  chapter  will  provide  an  overview  of  how  to  develop  complex  3-D  finite 
element  meshes  of  the  as-built  structures.  It  will  also  give  an  example  of 
how  to  mesh  around  a  crack  tip  to  determine  the  J-integral.  It  concludes 
by  showing  the  calculation  of  the  stress  intensity  factor  and  comparing  it 
against  the  material  toughness  (Kic). 

4.1  Meshing  techniques 

Historically  the  finite  element  mesh  generation  has  been  one  of  the  most 
time-consuming  issues  during  a  finite  element  analysis.  It  is  even  more 
complex  if  a  3-D  model  is  to  be  analyzed.  In  the  last  couple  of  years  new 
meshing  techniques  have  been  developed  to  construct  meshes  of  as-built 
structures  easily.  These  techniques  allow  the  modeler  to  include  the  actual 
details  of  the  structures  without  the  need  for  simplification  and  assump¬ 
tions.  Figure  16  shows  the  3-D  mesh  of  an  as-built  miter  gate.  The  inset 
provides  a  magnification  of  the  bottom  right  corner  of  the  gate.  It  becomes 
clear  that  advanced  meshing  techniques  are  necessary  to  generate  all  the 
details  shown  in  Figure  16  that  will  result  in  more  accurate  results. 

The  three  steps  to  generate  the  more  detailed  finite  element  mesh  are  as 
follows: 

•  A  3-D  CAD  drawing  of  the  as-built  structure  is  developed. 

•  Once  the  drawing  has  been  generated,  an  Initial  Graphics  Exchange 
Specifications  (IGES)  file  with  trimmed  surfaces  is  created.  Most  of  the 
commercial  CAD  packages  have  the  capability  to  generate  IGES  files. 

•  The  IGES  file  is  imported  into  the  finite  element  application  and  the 
mesh  is  generated. 

With  this  process  all  the  structural  components  are  easily  modeled  in  the 
CAD  package  instead  of  in  the  finite  element  mesh  application.  This 
utilizes  the  strengths  of  both  types  of  software;  the  CAD  program  to  create 
the  complex  geometry  and  the  finite  element  program  to  generate  the 
analytical  mesh  automatically,  making  the  meshing  inexpensive. 
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Figure  16.  3-D  Finite  element  mesh  of  as-built  miter  gate. 


For  this  example  the  finite  element  mesh  has  364,736  nodes  and  123,193 
shell  elements.  Each  shell  element  has  a  quadratic  formulation  with  five 
degrees  of  freedom  per  node.  The  total  number  of  degrees  of  freedom  for 
this  model  is  1.823  million. 

4.2  Load  and  displacement  boundary  condition 

For  this  example  the  load  boundary  condition  was  a  hydrostatic  head  of 
61.5  ft,  and  the  displacement  boundary  conditions  restrained  the  miter 
and  quoin  ends  to  move  along  the  three  principal  axes.  The  movement 
along  the  vertical  axis  at  the  pintle  was  also  restricted,  while  the  movement 
along  the  longitudinal  and  transverse  direction  (1  and  2)  was  restrained  at 
the  gudgeon.  Table  3  summarizes  the  boundary  conditions  used  for  this 
analysis. 


Table  3.  Boundary  conditions  used  for  the  fracture  mechanics  analysis. 


Gate  Component 

U1 

U2 

U3 

UR1 

UR2 

UR3 

CSYS 

Gugeon  pin 

0 

0 

— 

— 

— 

— 

Global 

Pintle 

— 

— 

0 

— 

— 

— 

Global 

Quoin  block 

0 

0 

0 

— 

— 

— 

Local 

Miter  block 

0 

0 

0 

— 

— 

— 

Local 
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4.3  Stress  intensity  factor 

The  downstream  flange  of  a  vertical  diaphragm  is  shown  in  Figure  17a,  and 
a  horizontal  girder  downstream  flange  is  shown  in  Figure  17b.  These 
flanges  meshes  contained  artificial  diagonal  cracks.  The  meshing  required 
to  calculate  the  J-integral  is  as  follows: 

•  The  crack  tip  must  be  meshed  with  a  ring  of  quarter-point  triangular 
elements. 

•  All  the  other  elements  outside  the  ring  can  be  quadratic  elements. 


Figure  17.  Finite  element  mesh  of  cracked  upstream  flange. 

Figure  18a  shows  the  stress  contour  of  the  upstream  vertical  diaphragm 
flange.  The  principal  stress  at  the  crack  tip  is  1,215  psi,  and  the  J-integral  is 
6.91  psi*in.  Using  Equation  2.15  and  solving  for  Ki  produced  a  stress 
intensity  factor  of  14.397 psi\ftn  ,  which  is  much  lower  than  the  fracture 
toughness  of  the  steel  (X/c=  50,000 psiyfin. ). 
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Figure  18b  shows  the  stress  contour  of  the  upstream  flange  of  the  hori¬ 
zontal  girder.  The  principal  stress  at  the  crack  tip  is  3,990  psi  and  the 
J-integral  is  24.66  psi*in.  Using  Equation  2.15  and  solving  for  Ki  produced 
a  stress  intensity  factor  of  27,202  psi-Jin. ,  which  is  much  lower  than  the 
fracture  toughness  of  the  steel  ( Kic=  50,000  psi-Jin. ). 

Since  the  calculated  stress  intensity  factors  are  less  than  the  fracture 
toughness  of  the  steel,  the  cracks  on  the  upstream  flanges  of  the  horizontal 
girder  and  vertical  diaphragm  will  not  propagate.  Whether  a  crack  will  or 
will  not  propagate  can,  therefore,  be  predicted  accurately.  Furthermore, 
ways  of  repairing  the  crack  regions  numerically  can  also  be  studied. 
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5  Numerical  Investigation  of  Diagonals  Pre 
Stress 


This  chapter  presents  a  comprehensive  3-D  finite  element  analysis  con¬ 
ducted  to  assess  the  diagonal  prestressing  levels.  This  analysis  will  provide 
the  initiation  for  a  new  design  approach.  The  new  techniques  will  allow  the 
prediction  of  optimum  prestressing  levels  without  compromising  the 
fatigue  capacity. 

In  recent  years  there  have  been  a  number  of  failures  of  diagonals  on  miter 
gates  (Figure  1).  These  appear  to  be  fatigue-induced  failures  driven  by  the 
connection  details  at  the  diagonal  ends,  but  some  evidence  also  suggests 
that  the  current  design  guidance  results  in  a  much  larger  prestress  in  miter 
gate  diagonals  than  may  be  required.  The  high  prestress  theoretically 
increases  the  fatigue  of  these  members  and  results  in  increased  stress  in 
the  miter  gate. 

A  series  of  parametric  analysis  were  carried  out  using  the  3-D  finite 
element  models  shown  in  Figure  19.  The  model  consisted  of  3-D  shell 
elements  with  6  degrees  of  freedom  (Figure  19). 


Downstream  Face 


Figure  19.  3-D  finite  element  model. 


The  following  parametric  cases  were  performed: 

1.  Gravity  with  no  diagonals:  gravity  analysis  of  the  miter  gate  without 
diagonals. 

2.  Gravity  with  positive  diagonal:  gravity  analysis  of  the  miter  gate  including 
only  the  positive  diagonal.  A  variation  of  the  diagonal  prestressing  was  also 
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analyzed  in  order  to  obtain  displacement  versus  diagonal  prestressing 
curves. 

3.  Gravity  with  positive  and  negative  diagonals:  gravity  analysis  of  the  miter 
gate  including  the  positive  and  negative  diagonals.  Variation  in  the 
diagonal  prestressing  was  also  carried  out. 

4.  Gravity  with  positive  and  negative  diagonals,  hydrostatic  resistance 
pressure,  and  movement:  analysis  of  the  miter  gate  including  the  positive 
and  negative  diagonals,  the  hydrostatic  resistance  pressure,  and  the 
movement  of  the  gate  were  evaluated.  Diagonal  prestressing,  hydrostatic 
resistance  pressure  loads,  and  movement  of  the  gate  with  different 
diagonals  prestressing  were  also  conducted. 

5.  Quoin  block  analysis:  gate  in  closed  position  with  upper  and  lower  pools.  A 
typical  combination  of  diagonal  prestressing  was  used  in  the  analysis.  The 
prestressing  combination  was  20  ksi  in  the  positive  diagonal  and  25  ksi  in 
the  negative  diagonal.  These  analyses  account  for  deterioration  in  the 
quoin  block  area. 

5.1  Diagonal  prestressing 

A  method  to  produce  the  same  effect  of  the  diagonal  prestressing  in  the 
finite  element  model  when  using  plate  elements  was  evaluated.  In  order  to 
represent  the  pre-stress  in  the  diagonals  using  the  finite  element  code,  a 
thermal  gradient  (AT)  of  the  diagonals  was  used.  To  determine  the  sign  of 
the  change  in  temperature,  a  shell  element  was  created  to  model  a  bar.  The 
bar  was  assumed  to  be  fixed  at  its  ends.  To  compare  the  results  obtained 
with  the  model,  a  theoretical  solution  of  the  problem  was  used.  The 
theoretical  formula  for  a  fully  restrained  element  is: 

o  =  -Ecl(AT)  (5.1) 

Where  E  is  Young’s  modulus  of  elasticity,  a  measures  the  percentage 
change  in  the  length  per  degree  temperature  (Equation  5.2).  In  equation 
5.2,  L0  is  the  original  length  and  AL  is  the  change  in  length  after 
temperature  gradient  ( AT)  is  applied.  For  steel  a  =  6.5  e-6  /  °F. 


f  AL 


AT 


(5.2) 
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Figure  20  presents  a  free  body  diagram  of  the  pre-stressed  element.  In 
order  to  produce  a  tensile  stress  in  the  member,  a  negative  AT  is  required. 
The  negative  AT  produces  contraction  (shortening)  in  the  element,  this 
contraction  in  a  restrained  element  produce  stresses  similar  to  the  ones 
produced  by  the  prestressing  force  in  the  miter  gate. 


Figure  20.  Free  body  diagram  of  the  pre¬ 
stressed  element. 


A  2-  by  2-  in.  steel  section  with  E  =  30  x  106  psi  and  different  temperature 
gradients  was  used.  The  theoretical  and  finite  element  results  with 
different  temperature  gradients  are  compared  in  Table  4. 

Table  4.  Validation  of  the  analysis  using  different  values  of  AT. 


A  T 

0,  ksi 

Theoretical 

Finite  Element 

Model 

-20 

3.9 

3.9 

-40 

7.8 

7.8 

-60 

11.7 

11.7 

-80 

15.6 

15.6 

5.2  Automatic  stabilization  of  static  problems 

In  finite  elements  the  analyses  that  include  nonlinearities  can  experience 
instabilities.  The  instabilities  may  be  of  a  geometrical  nature  such  as 
buckling  or  a  material  nature  such  as  material  softening.  If  the  instabil¬ 
ities  manifest  themselves  in  a  global  load-displacement  response  with  a 
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negative  stiffness,  the  problem  can  be  treated  as  a  buckling  or  collapse 
problem.  But  if  the  instability  is  localized,  a  local  transfer  of  strain  energy 
from  one  part  of  the  model  to  neighboring  parts  occurs;  and  global 
solution  methods  may  not  work  properly.  These  problems  have  to  be 
solved  either  dynamically  or  with  the  aid  of  artificial  damping.  The 
commercial  code  used  provides  an  automatic  mechanism  for  stabilizing 
unstable  quasi-static  problems  through  the  automatic  addition  of  volume 
proportional  damping  to  the  model.  Viscous  forces  with  the  form 
Fv  =  cMv  where  c  is  a  damping  factor,  M  =  artificial  mass  matrix 

calculated  with  unity  density,  and  v  =  An/ At  is  the  vector  of  nodal  velocity, 
are  added  to  the  global  equilibrium  equations  P  — I  — Fv=  0  where  P  is  the 

external  force,  and  I  is  the  internal  force. 

Because  of  the  complexity  of  the  model,  localized  nonlinearities  can  occur 
without  affecting  the  global  stability  of  the  structure.  For  that  reason  this 
automatic  stabilization  mechanism  was  activated  in  the  parametric 
analysis.  The  value  used  for  the  damping  factor  was  the  default  value  of 
0.0002. 

5.3  Miter  gate  operation 

During  gate  operation  the  boundary  conditions  are  located  in  the  pintle 
and  gudgeon  pin.  These  boundary  conditions  do  not  restrict  the  gate  from 
having  a  rigid  body  rotation  around  the  axis  connecting  the  pintle  and 
gudgeon  pin.  To  induce  a  rotation  in  the  gate  three  alternatives  were 
evaluated:  (l)  applied  rotation  to  the  pintle,  (2)  applied  angular  velocity  to 
the  pintle,  and  (3)  a  concentrated  force  in  the  location  of  the  operation 
point.  The  load  conditions  for  each  alternative  were  the  self-weight  of  the 
gate,  prestressing  in  the  diagonals,  and  the  water  resistance  pressure 
obtained  from  code  recommendations,  30  psf  (Headquarters,  U.S.  Army 
Corps  of  Engineers  1994).  The  analysis  reveals  that  the  three  methods 
induce  rotation  in  the  gate. 

The  applied  rotation  in  the  pintle  and  the  angular  velocity  cases  generated 
required  rotary  motion  in  the  gate.  However,  the  relative  displacements 
between  the  top  and  bottom  of  the  gate  induced  by  the  operating  force 
were  not  achieved.  The  third  alternative  (concentrated  force)  produced  a 
more  realistic  behavior  of  the  gate  when  in  operation;  however,  because 
the  rotation  cannot  be  specified,  the  total  rotation  depends  on  the  step 
time.  To  reproduce  the  required  rotation  during  the  gate  operation,  the 
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rigid  body  rotation  had  to  be  restrained.  To  restrain  the  rotation  of  the 
gate,  contact  surfaces  were  created.  The  contact  surfaces  act  as  a 
restriction  in  the  model,  stopping  the  gate  at  the  desired  locations.  In 
order  to  verify  the  effectiveness  of  the  method,  a  simplified  model  was 
created  with  the  same  boundary  condition  of  the  real  structure.  In  this 
model  a  60-  by  40-ft  plate  was  used  to  represent  the  gate  and  two  plates 
were  used  to  act  as  restricting  walls.  Figure  21  presents  a  view  of  the  model 
rotating  between  the  restricting  elements. 


Figure  21.  Simplified  model  master  and  slave  selection. 


To  make  the  gate  stop  at  the  desired  location  a  surface-to-surface  contact 
interaction  between  the  surfaces  was  specified.  Surface-to-surface  contact 
interactions  describe  contact  between  two  deformable  surfaces  or  between 
a  deformable  surface  and  a  rigid  surface.  The  finite  element  program 
defines  contact  between  two  bodies  in  terms  of  two  surfaces  that  may 
interact;  these  surfaces  are  called  a  “contact  pair.”  The  contact  surfaces 
are  defined  as  master  and  slave  surfaces.  The  surface  of  the  stiffer  body 
should  be  chosen  as  the  master  surface,  and  the  slave  surface  should 
always  be  attached  to  deformable  bodies. 

The  master  surface  normal  always  has  to  be  in  the  direction  of  the  contact 
surface.  In  this  case  the  master  surfaces  are  the  two  shell  elements  located 
outside  of  the  model  and  the  slave  surface  is  the  gate  (Figure  21). 
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After  the  selection  of  the  master  and  slave  surfaces,  a  contact  property 
needs  to  be  defined.  The  relationship  of  default  contact  pressure  over 
closure  used  by  the  code  is  the  “hard”  contact  model.  This  model  assumes 
that  the  surfaces  transmit  no  contact  pressure  unless  the  nodes  of  the  slave 
surface  contact  the  master  surface.  There  is  no  limit  to  the  magnitude  of 
contact  pressure  that  can  be  transmitted  when  the  surfaces  are  in  contact. 
That  means  that  the  slave  surface  cannot  penetrate  the  master  surface. 
Because  the  surfaces  make  contact  in  a  localized  area,  the  mechanical 
property  used  was  “normal  behavior.”  The  small  sliding  formulation  was 
specified  in  the  analysis  because  the  master  and  slave  surfaces  are  in 
contact  in  a  specific  area.  The  formulation  assumes  that  the  surfaces  may 
undergo  arbitrarily  large  rotations  but  that  a  slave  node  will  interact  with 
the  same  local  area  of  the  master  surface  throughout  the  analysis.  The 
formulation  used  in  this  simplified  model  was  implemented  with  success 
in  the  miter  gate  3-D  model  to  produce  the  desired  behavior. 

5.4  Augmented  Lagrangian  method 

To  permit  some  penetration  of  the  slave  surface  into  the  master  surface, 
the  augmented  Lagrangian  method  was  used.  If  a  slave  node  penetrates 
the  master  surface  by  more  than  0.1  percent  of  the  characteristic  interface 
length,  the  contact  pressure  is  augmented  and  another  series  of  iterations 
is  executed  until  convergence  is  once  again  achieved.  Only  when  the  pene¬ 
tration  tolerance  requirement  is  satisfied  will  the  solution  be  accepted. 

This  method  may  require  additional  iterations  in  some  cases;  however, 
this  approach  is  sometimes  useful  for  numerical  reasons  because  it  can 
make  the  resolution  of  contact  conditions  easier. 

5.5  Parametric  analysis  results 

This  section  presents  the  parametric  analysis  results  for  the  different 
analyzed  cases.  The  vertical  and  lateral  displacements  were  measured  in 
the  miter  gate  end  in  the  top  and  bottom  locations  for  all  the  cases.  The 
results  presented  herein  summarize  the  relevant  findings  in  the  different 
cases  performed. 

5.5.1  Case  1:  Gravity  with  no  diagonals 

This  analysis  was  carried  out  to  study  the  behavior  of  the  miter  gate  when 
subjected  to  the  self  weight  (gravity  effects)  only.  Figure  22  shows  a 
schematic  drawing  of  this  case.  The  boundary  conditions  used  are 
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presented  in  Table  5.  Because  the  pintle  and  the  gudgeon  pin  restrict  one 
end  of  the  miter  gate,  the  results  demonstrate  that  the  gate  behaves  like  a 
cantilever  beam,  producing  tension  in  the  gudgeon  pin  zone  and  com¬ 
pression  at  the  pintle  zone.  Figure  23  shows  the  vertical  displacement  of 
the  miter  end.  During  the  gravity  analysis  the  results  revealed  a  twist  in 
the  miter  gate  end.  The  quoin  end  of  a  miter  gate  is  held  vertical  by  the 
pintle  and  gudgeon  pin,  leaving  the  miter  end  free  to  twist  out  of  vertical 
alignment.  This  twist  is  produced  because  the  center  of  gravity  and  the 
shear  center  are  not  in  the  same  vertical  plane.  This  difference  in  its  vert¬ 
ical  planes  produces  an  eccentricity  that  produces  a  load-displacement 
(P-Delta)  effect  in  the  structure.  Figure  24  presents  the  observed  behavior 
with  more  details.  The  lateral  displacement  for  this  case  is  presented  in 
Figure  25. 


Table  5.  Boundary  condition  for  the  analyzed  structure. 


Gate 

Component 

U1 

U2 

U3 

UR1 

UR2 

UR3 

CSYS 

Gudgeon  pin 

0 

0 

— 

— 

— 

— 

Global 

Pintle 

0 

0 

0 

0 

0 

— 

Global 

Quoin  block 

— 

— 

0 

— 

— 

— 

Local 

Miter  block 

— 

— 

— 

— 

— 

— 

— 
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Figure  23.  Vertical  displacement. 


Figure  24.  Lateral  observed  behavior. 
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Gravity  With  No  Diagonals  Analysis 
Lateral  Displacement  [U2G] 

Downstream  Displacement  [ft]  Upstream 


Figure  25.  Lateral  displacement. 

5.5.2  Case  2:  Gravity  with  positive  diagonals 

To  study  the  effects  of  the  positive  diagonal  in  the  structure,  a  series  of 
analyses  were  conducted  including  the  gravity  loads  with  different  positive 
diagonal  prestressing.  Figure  26  presents  a  schematic  drawing  of  the 
structure  used  in  this  case.  The  boundary  conditions  used  are  presented  in 
Table  6.  The  data  presented  in  Figures  27  and  28  and  Table  8  correspond 
to  top  (Point  A)  and  bottom  (Point  B)  nodes  at  the  miter  end.  Table  7 
shows  the  different  diagonal  pre-stressing  combinations  used  for  the  posi¬ 
tive  diagonal.  Table  8  summarizes  the  spatial  displacements  obtained  for 
the  different  diagonal  prestressing  used  in  the  analysis. 

Vertical  displacement: 

The  analysis  demonstrated  that  the  structure  experienced  a  negative  vert¬ 
ical  displacement  (downward  displacement)  because  of  the  self  weight  of 
the  miter  gate.  To  reduce  the  excessive  deflection  during  operation,  diag¬ 
onals  were  used  in  the  structure.  The  positive  diagonal  pulled  up  the  miter 
gate;  thus,  the  greater  the  prestressing  in  the  positive  diagonals,  the 
smaller  the  vertical  deflection  became. 
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Table  6.  Boundary  conditions  for  the  analyzed  structure. 


Gate  Component 

U1 

U2 

U3 

UR1 

UR2 

UR3 

CSYS 

Gudgeon  pin 

0 

0 

— 

— 

— 

— 

Global 

Pintle 

0 

0 

0 

0 

0 

— 

Global 

Quoin 

— 

— 

0 

— 

— 

— 

Local 

Figure  27.  Variation  in  vertical  displacement  for  the  different  diagonal 

prestressing. 
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Pos.  Diag.  Pre-stressing  vs  Displacement 
U3  (Top) 


£ 

&J0 


f/3 

f/5 

£ 

Ph 


Displacement  [ft] 


Figure  28.  Top  vertical  displacement  variation. 


Table  7.  Positive  diagonal  prestressing. 


Case 

Diagonal  Prestressing,  ksi  (Positive) 

1 

0 

2 

20 

3 

40 

4 

60 

5 

80 

Table  8.  Spatial  displacement  summary  for  the  top  and  bottom  nodes  of  the  miter  gate. 


Displacement,  ksi 

Node 

U1  [ft] 

U2  [ft] 

U3  [ft] 

0 

Top  (A) 

0.00064 

-0.09772 

-0.13240 

Bottom  (B) 

-0.00457 

0.08347 

-0.13225 

20 

Top  (A) 

0.00073 

0.04640 

-0.09660 

Bottom  (B) 

-0.00678 

-0.04169 

-0.09630 

40 

Top  (A) 

0.00051 

0.19160 

-0.06210 

Bottom  (B) 

-0.00927 

-0.16670 

-0.06038 

60 

Top  (A) 

-0.00006 

0.33120 

-0.02802 

Bottom  (B) 

-0.01207 

-0.29020 

-0.02402 

80 

Top  (A) 

-0.00089 

0.48110 

0.00609 

Bottom  (B) 

-0.01496 

-0.41680 

0.01371 
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The  vertical  displacement  became  positive  for  a  certain  prestress  value. 
The  prestress  required  to  produce  zero  vertical  displacement  was  greater 
than  the  allowable  stress  of  the  material,  which  was  60  ksi. 

Figure  27  presents  the  variation  in  vertical  displacement  for  the  different 
diagonal  prestressing.  Figure  28  presents  in  graphical  form  the  vertical 
displacement  variation  for  the  top  node  in  this  analysis.  Using  these 
graphs  the  diagonal  prestressing  in  the  positive  diagonal  required  to 
produce  zero  vertical  displacement  can  be  obtained,  about  70  ksi. 

Lateral  displacement. 

In  this  analysis  the  self  weight  of  the  miter  gate  and  the  positive  diagonal 
prestressing  caused  the  gate  to  twist.  In  this  step  the  positive  diagonal 
contribution  to  the  structure  in  the  lateral  displacements  was  analyzed. 

The  analysis  reveals  an  upstream  movement  of  the  top  portion  of  the  miter 
end  (Point  A)  with  an  increase  in  the  positive  diagonal  prestressing.  The 
greater  the  prestressing,  the  bigger  the  upstream  twist  became.  Figure  29 
presents  the  observed  behavior  for  the  different  diagonal  prestressing.  The 
gravity  produced  a  downstream  movement  of  the  miter  end,  represented 
in  the  curve  in  which  the  diagonal  prestressing  was  o  ksi,  and  the  positive 
diagonal  prestressing  produced  an  upstream  movement  of  the  top  portion 
of  the  gate  end  (Point  A)  when  the  prestressing  increased  as  observed  in 
the  plot.  Figure  30  shows  that  a  diagonal  prestressing  of  14  ksi  will 
eliminate  any  relative  displacement  between  Point  A  and  Point  B. 

5.5.3  Case  3:  Gravity  with  positive  and  negative  diagonals 

Case  3  includes  the  gravity  loads  and  the  positive  and  negative  diagonals. 
Different  diagonal  prestressing  combinations  were  used  in  this  case.  The 
selected  prestressing  values  were  o,  5, 10, 15,  20,  25,  and  30  ksi.  Figure  31 
presents  a  schematic  view  of  the  structure  for  this  case.  The  boundary  con¬ 
ditions  used  are  offered  in  Table  9.  Figures  32-37  and  Tables  10  and  11 
present  a  summary  of  the  spatial  displacements  at  the  top  and  bottom  of 
the  miter  end. 
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Gravity  +  Pos.  Diagonal  Pre-stressing  Analysis 
Lateral  Displacement  [U2G] 

Downstream  Displacement  [ft]  Upstream 


"0  ksi 


20  Ksi 


40  ksi 


60  ksi  )K  80  ksi 


Figure  29.  Lateral  behavior  of  the  miter  gate. 


Pos.  Diag.  Pre-stressing  vs.  Relative  Lateral 
Displacement 


Relative  Lateral  Displacement  [ft] 


Figure  30.  Relative  lateral  displacement. 
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Table  9.  Boundary  condition  for  the  analyzed  structure. 


Gate 

Component 

U1 

U2 

U3 

UR1 

UR2 

UR3 

CSYS 

Gudgeon  pin 

0 

0 

— 

— 

— 

— 

Global 

Pintle 

0 

0 

0 

0 

0 

— 

Global 

Quoin  block 

— 

— 

0 

— 

— 

— 

Local 

Miter  block 

— 

— 

— 

— 

— 

— 

Vertical  displacement: 

The  analysis  confirmed  that  the  structure  had  a  negative  vertical  displace¬ 
ment  caused  by  the  self  weight  (gravity  loads)  of  the  miter  gate.  Diagonals 
helped  reduce  the  excessive  deflection  during  operation.  The  greater  the 
prestressing  in  the  positive  diagonals,  the  smaller  the  vertical  deflection 
became.  The  negative  diagonals  had  the  inverse  effect:  the  greater  the 
prestressing  in  the  negative  diagonals,  the  greater  the  vertical  displace¬ 
ment  became.  Table  10  shows  the  vertical  displacement  results  for  the 
different  diagonal  prestressing  combinations.  Results  showed  that  the 
minimum  vertical  displacement  of  0.0007  ft  occurred  with  10  ksi  of 
prestressing  in  the  negative  diagonal  and  20  ksi  in  the  positive  diagonal. 
Figure  32  presents  the  variation  in  vertical  displacement  of  the  miter  gate 
end  for  the  different  combination  of  diagonal  prestressing.  Each  series 
corresponds  to  a  positive  diagonal  prestressing,  and  the  vertical  axis 
represents  the  negative  diagonal  prestressing.  The  horizontal  axis 
represents  the  vertical  displacement  of  the  miter  end. 
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Table  10.  Vertical  displacement  summary. 


Positive  Diagonal 

0 

5 

10 

15 

20 

25 

30 

Negative  Diagonal 

0 

-0.01200 

-0.00606 

-0.00092 

0.00342 

0.00775 

0.01209 

0.01644 

5 

-0.01302 

-0.00870 

-0.00440 

-0.00009 

0.00423 

0.00857 

0.01291 

10 

-0.01649 

-0.01223 

-0.00794 

-0.00363 

0.00070 

0.00503 

0.00938 

15 

-0.02005 

-0.01578 

-0.01148 

-0.00717 

-0.00285 

0.00149 

0.00584 

20 

-0.02362 

-0.01934 

-0.01505 

-0.01073 

-0.00640 

-0.00207 

0.00229 

25 

-0.02719 

-0.02291 

-0.01861 

-0.01429 

-0.00996 

-0.00562 

-0.00127 

30 

-0.03077 

-0.02648 

-0.02218 

-0.01786 

-0.01352 

-0.00919 

-0.00483 

Gravity  +  Neg.  and  Pos.  Diag.  Prestressing 
Vertical  Displacement  [U3G] 


Displacement  [ft] 


♦  0  ksi  ■  5  ksi  A  10  ksi  15  ksi  )l(  20  ksi  •  25  ksi  I  30  ksi 

Figure  32.  Vertical  displacement  variation  for  the  different  diagonal  prestressing. 

In  Figure  32  the  negative  diagonal  prestressing  required  to  obtain  zero 
vertical  displacement  is  the  intersection  between  the  lines  with  the  vertical 
axis.  For  example,  for  20  ksi  of  prestressing  in  the  positive  diagonal,  11  ksi 
of  prestressing  in  the  negative  diagonal  is  required  to  have  zero  vertical 
displacement  in  the  miter  end. 
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Lateral  displacement: 

The  positive  and  negative  diagonal  prestressing  contributions  in  the  lateral 
displacement  were  analyzed.  The  self  weight  of  the  miter  gate  produces  a 
downstream  movement  of  the  top  portion  of  the  miter  end  (Figure  25). 

The  positive  diagonal  prestressing  was  evaluated  in  the  previous  case,  and 
the  behavior  was  an  upstream  movement  of  the  top  portion  of  the  miter 
end  as  the  prestressing  increased,  while  the  bottom  portion  of  the  miter 
gate  moved  downstream.  Table  11  summarizes  the  lateral  displacements 
obtained  for  the  different  positive  and  negative  diagonal  prestressing 
combinations.  Figure  33  presents  a  typical  plot  of  the  observed  lateral 
behavior  of  the  miter  end  with  a  positive  diagonal  prestressing  of  15  ksi, 
and  varying  negative  diagonal  prestressing.  Using  this  plot,  the  negative 
diagonal  effect  in  the  lateral  displacement  of  the  structure  can  be  obtained. 
The  plot  shows  an  initial  location  of  the  top  portion  of  the  miter  end  in  the 
upstream  position  caused  by  the  positive  diagonal  prestressing  of  15  ksi. 


Table  11.  Spatial  displacement  summary  for  the  top  and  bottom  of  the  miter  gate. 


Top 

Positive  Diagonal 

0 

5 

10 

15 

20 

25 

30 

Negative  Diagonal 

0 

-0.03499 

-0.01057 

0.00966 

0.02556 

0.04140 

0.05725 

0.07313 

5 

-0.03985 

-0.02443 

-0.00872 

0.00703 

0.02284 

0.03866 

0.05450 

10 

-0.05863 

-0.04303 

-0.02735 

-0.01161 

0.00417 

0.01999 

0.03583 

15 

-0.07733 

-0.06173 

-0.04605 

-0.03031 

-0.01452 

0.00130 

0.01712 

20 

-0.09610 

-0.08048 

-0.06479 

-0.04905 

-0.03329 

-0.01743 

-0.00159 

25 

-0.11490 

-0.09927 

-0.08355 

-0.06780 

-0.05203 

-0.03621 

-0.02035 

30 

-0.13373 

-0.11806 

-0.10234 

-0.08660 

-0.07081 

-0.05497 

-0.03912 

Bottom 

Positive  Diagonal 

0 

5 

10 

15 

20 

25 

30 

Negative  Diagonal 

0 

0.02266 

-0.00208 

-0.02472 

-0.04518 

-0.06568 

-0.08621 

-0.10676 

5 

0.02547 

0.00507 

-0.01527 

-0.03569 

-0.05617 

-0.07669 

-0.09725 

10 

0.03488 

0.01464 

-0.00571 

-0.02614 

-0.04663 

-0.06715 

-0.08772 

15 

0.04453 

0.02425 

0.00387 

-0.01656 

-0.03705 

-0.05759 

-0.07817 

20 

0.05420 

0.03388 

0.01349 

-0.00696 

-0.02747 

-0.048 

-0.06859 

25 

0.06388 

0.04353 

0.02312 

0.00265 

-0.01786 

-0.03842 

-0.05901 

30 

0.07355 

0.05319 

0.03276 

0.01227 

-0.00826 

-0.0288 

-0.04942 
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An  increase  in  the  negative  diagonal  prestressing  produces  a  downstream 
movement  of  the  top  portion  of  the  miter  end  and  an  upstream  movement 
of  the  bottom  portion  of  the  miter  end.  This  effect  is  produced  by  the  pre¬ 
stressing  in  the  positive  and  negative  diagonals.  The  greater  the  pre¬ 
stressing,  the  bigger  the  twist  becomes  in  that  direction. 

Figure  34  presents  top  (TD)  and  bottom  (BD)  lateral  displacements  for  all 
49  analyses  conducted.  The  plot  shows  the  miter  end  behavior  varying  the 
negative  diagonal  prestressing  as  the  y-axis  versus  the  lateral  displace¬ 
ment.  Each  series  corresponds  to  a  constant  positive  diagonal  pre¬ 
stressing.  The  series  with  negative  slope  represents  the  displacements  in 
the  top  of  the  miter  end,  and  the  series  with  positive  slope  represents  the 
displacements  in  the  bottom  of  the  miter  end. 
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Gravity  +  Neg.  and  Pos.  Diag.  Prestressing 
Lateral  Displacement  [U2G] 


TD  TD  TD  TD  TD  BD  BD  BD  BD  BD 

Oksi  5  ksi  lOksi  15  ksi  20ksi  20ksi  15  ksi  lOksi  5  ksi  Oksi 


Downstream  Displacement  [ft]  Upstream 


"0  ksi  (Top) 
10  ksi  (Bot) 


5  ksi  (Top)  -±- 10  ksi  (Top)  15  ksi  (Top)  20  ksi  (Top)  — ' ■—  0  ksi  (Bot)  — ■—  5  ksi  (Bot) 

15  ksi  (Bot)  — —  20  ksi  (Bot)  25  ksi  (Top)  30  ksi  (Top)  25  ksi  (Bot)  30  ksi  (Bot) 


Figure  34.  Positive  and  negative  diagonals  relation. 


From  Figure  34  the  required  positive  and  negative  diagonal  prestressing 
necessary  to  obtain  zero  relative  displacement  between  the  top  and  bottom 
ends  is  given  by  the  point  of  intersection  of  series  representing  top  and 
bottom  displacements  for  the  same  positive  diagonal  prestressing.  This 
point  of  intersection  is  associated  with  a  downstream  displacement.  The 
downstream  displacement  occurs  because  the  diagonals  are  located  in  the 
downstream  side  of  the  gate,  and  an  increase  in  the  prestressing  of  the 
diagonals  produces  a  small  rotation  of  the  gate  in  the  downstream  direc¬ 
tion.  This  rotation  is  less  than  1  in.  in  all  the  prestressing  combinations.  An 
approximation  of  the  top  and  bottom  displacements  can  also  be  obtained 
for  the  different  combinations  of  diagonal  prestressing.  Figure  35  shows 
that  for  a  negative  diagonal  prestressing  of  15  ksi  and  a  positive  diagonal 
prestressing  of  0.0  ksi,  the  top  of  the  gate  will  displaced  0.077  ft  down¬ 
stream  and  the  bottom  will  move  0.044  ft  upstream.  Similar  evaluations 
can  be  accomplished  for  different  combinations  of  diagonal  prestressing. 

Figures  36  and  37  present  the  relative  lateral  displacements  in  which  a 
variation  of  the  negative  diagonal  is  applied  with  a  constant  positive 
diagonal  prestressing  and  when  a  variation  of  the  positive  diagonal  is 
applied  with  a  constant  negative  diagonal  prestressing,  respectively. 
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Gravity  +  Neg.  and  Pos.  Diag.  Prestressing 
Lateral  Displacement  [U2G] 


TD  TD  BD  BD 

Oksi  15ksi  15ksi  0  ksi 


Downstream  Displacement  [ft]  Upstream 


"0  ksi  (Top) 


15  ksi  (Top)  1  0  ksi  (Bot) 


15  ksi  (Bot) 


Figure  35.  Detailing  of  the  positive  and  negative  diagonal  relation. 


Gravity  +  Neg.  and  Pos.  Diag.  Prestressing 
Relative  Lateral  Displacement  [U2G] 


Downstream  Displacement  [ft]  Upstream 


0  ksi  ■  5  ksi  10  ksi  15  ksi  )K  20  ksi  •  25  ksi  I  30  ks: 


Figure  36.  Relative  lateral  displacement  (negative  diagonals  vary). 
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Gravity  +  Neg.  and  Pos.  Diag.  Prestressing 
Relative  Lateral  Displacement  [U2G] 


Downstream  Displacement  [ft]  Upstream 


♦  0  ksi  ■  5  ksi  A  10  ksi 

15  ksi  X  20  ksi  •  25  ksi  1  30  ksi 

Figure  37.  Relative  lateral  displacement  (positive  diagonals  vary). 


These  plots  can  be  used  to  determine  the  required  combination  of  diagonal 
prestressing  to  produce  zero  relative  lateral  displacement  in  the  miter  end. 
The  series  represents  the  constant  diagonal  prestressing.  Using  Figure  36 
the  series  represents  the  positive  diagonal  prestressing,  and  the  vertical 
axis  represents  the  negative  diagonal  prestressing.  When  a  series  inter¬ 
sects  the  vertical  axis,  the  point  of  intersection  represents  the  negative 
diagonal  prestressing  required  to  produce  zero  relative  displacement  in 
the  miter  end.  For  example  from  Figure  36  for  15  ksi  of  prestressing,  in  the 
positive  diagonal  12.5  ksi  of  prestressing  is  required  in  the  negative 
diagonal  to  produce  zero  relative  displacement  in  the  gate. 

5.5.4  Case  4:  Gravity  with  positive,  negative  diagonals,  movement,  and 
water  resistance 

This  parametric  case  was  directed  to  study  the  behavior  of  the  gate  when  it 
is  opening  and  closing.  A  total  of  16  analyses  were  performed.  In  Case  3 
the  behavior  of  the  gate  miter  end  was  evaluated  for  gravity  and  different 
diagonal  prestressing.  In  addition  to  the  previous  case,  the  effects  of  the 
water  resistance  pressures  and  the  gate  movement  were  evaluated.  The 
analysis  was  divided  into  two  steps,  opening  the  gate  and  closing  the  gate. 
Because  the  water  levels  upstream  and  downstream  of  the  miter  gate  when 
the  gate  is  in  operation  are  the  same,  the  analysis  considered  the  resis¬ 
tance  pressure  produced  by  the  water  only  when  the  gate  is  opening  and 
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closing.  EM  1110-2-2703  (Headquarters,  U.S.  Army  Corps  of  Engineers, 
1994)  and  EM  1110-2-2105  (Headquarters,  U.S.  Army  Corps  of  Engineers, 
1993)  suggest  water  resistance  pressures  should  be  either  30  psf  or  45  psf 
depending  on  such  gate  geometrical  properties  as  height  and  length.  For 
gates  similar  to  the  one  analyzed  in  this  study  the  code  suggests  a  value  for 
water  resistance  pressure  of  30  psf.  The  water  resistance  pressure  alone 
produces  downstream  movement  of  the  bottom  portion  of  the  gate  miter 
end  with  respect  to  its  neutral  axis;  the  top  portion  moves  in  the  upstream 
direction.  To  induce  movement  in  the  gate,  a  concentrated  force  was 
applied  at  the  top  of  the  gate  in  the  operation  point  of  the  gate.  In  order  to 
simulate  the  gate  opening  and  closing,  the  concentrated  force  applied  was 
represented  in  two  steps,  one  with  positive  magnitude  to  open  the  gate  and 
the  other  with  negative  magnitude  to  close  the  gate.  Because  of  the  large 
displacements  produced  for  the  gate  movement,  the  concentrated  force 
was  applied  following  the  rotation  of  the  applied  nodes.  A  concentrated 
force  of  250  to  300  kips  is  typically  applied  to  open  and  close  this  type  of 
gate.  In  this  study  the  load  magnitude  used  was  250  kips.  This  load 
applied  in  the  operation  point  produced  an  upstream  movement  of  the  top 
portion  of  the  gate  during  opening  and  a  downstream  movement  of  the  top 
portion  of  the  gate  during  closing.  Figure  38  shows  a  schematic  drawing  of 
the  gate  for  both  steps:  opening  and  closing.  Table  12  presents  the 
boundary  conditions  used  for  this  case. 

Vertical  displacement: 

The  previous  analysis  confirmed  that  the  self  weight  of  the  gate  and  the 
prestressing  of  the  different  diagonals  produce  changes  in  the  vertical 
displacements  of  the  miter  end.  In  this  case  the  vertical  behavior  of  the 
miter  end  was  evaluated  for  effects  of  the  water  resistance  and  opening 
and  closing  the  gate. 

The  results  reveal  an  upward  movement  of  the  gate  miter  end  during 
opening  and  a  downward  movement  of  the  gate  end  during  closing.  These 
values  are  less  than  0.6  in.,  negligible  in  comparison  with  the  gate 
dimensions. 
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Table  12.  Boundary  condition  for  the  analyzed  structure. 


Gate 

Component 

U1 

U2 

U3 

UR1 

UR2 

UR3 

CSYS 

Gudgeon  pin 

0 

0 

— 

— 

— 

— 

Global 

Pintle 

0 

0 

0 

0 

0 

— 

Global 

Quoin  block 

Miter  block 

_ 

_ 

_ 

_ 

— 

— 

Local 
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Table  13  presents  the  vertical  displacement  results  for  the  different 
diagonal  prestressing  combinations  when  the  gate  was  opening  and 
closing.  Figure  39  presents  the  vertical  displacements  obtained  for  each 
diagonal  prestressing  combination. 


Table  13.  Vertical  displacement  results. 


Positive  Diagonal 

0 

10 

20 

30 

Opening 

Negative 

Diagonal 

0 

0.0205 

0.0221 

0.0341 

0.0471 

10 

-0.0003 

0.0105 

0.0234 

0.0365 

20 

-0.0132 

-0.0003 

0.0128 

0.0259 

30 

-0.0240 

-0.0110 

0.0021 

0.0153 

Positive  Diagonal 

0 

10 

20 

30 

Closing 

Negative 

Diagonal 

0 

-0.0596 

-0.0360 

-0.0124 

0.0028 

10 

-0.0585 

-0.0348 

-0.0209 

-0.0077 

20 

-0.0578 

-0.0446 

-0.0315 

-0.0184 

30 

-0.0683 

-0.0553 

-0.0422 

-0.0290 

Gravity  +  Neg.  Pos.  Diag.  Pre-stressing  +  Movement  + 
Water  Resistance 
Vertical  Displacement 


OX 


<L> 

U 


QJ 

U 

Pm 

ox 

& 


OX 

CJ 

£ 


Displacement  [ft] 


♦  0  ksi  Open  1 
)K  20  ksi  Open 


•  0  ksi  Close  ' 
20  ksi  Close  ' 


•  10  ksi  Open 

•  30  ksi  Open 


10  ksi  Close 
30  ksi  Close 


Figure  39.  Vertical  displacements. 
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Lateral  displacement: 

The  boundary  conditions  of  a  moving  gate  are  located  only  in  the  pintle 
and  gudgeon  pin.  Therefore,  the  gate  can  rotate  freely,  with  the  only 
restriction  being  the  water  resistance  pressure.  The  self  weight  of  the  miter 
gate,  the  positive  and  negative  diagonals  prestressing,  and  the  movement 
of  the  gate  cause  it  to  twist.  In  these  cases  the  contributions  of  the  positive 
and  negative  diagonals  when  the  gate  is  opening  and  closing  were  ana¬ 
lyzed.  The  previous  analysis  revealed  that  prestressing  the  positive  diag¬ 
onals  produces  a  twist  in  the  miter  gate,  causing  the  top  portion  of  the 
miter  end  to  move  upstream;  and  prestressing  the  negative  diagonals 
produces  a  downstream  movement  of  the  top  portion  of  the  miter  end. 
These  sets  of  analysis  were  conducted  to  asses  the  behavior  of  the  miter 
end  when  the  gate  is  moving  and  the  contribution  of  the  different 
combinations  of  diagonals  prestressing. 

The  incorporation  of  the  gate  movement  produces  an  upstream  movement 
of  the  top  portion  of  the  gate  when  opening  and  a  downstream  movement 
of  the  top  portion  of  the  gate  when  closing.  The  critical  lateral  spatial  dis¬ 
placements  are  summarized  in  Table  14  for  the  cases  when  the  gate  is 
opening  and  in  Table  15  for  the  cases  when  the  gate  is  closing.  These 
maximum  values  occur  immediately  after  the  operation  load  is  applied. 


Table  14.  Lateral  displacement  results. 


Top  End 

Positive  Diagonal 

0 

10 

20 

30 

Negative 

Diagonal 

0 

0.1081 

0.1164 

0.1719 

0.2330 

10 

0.0118 

0.0625 

0.1234 

0.1849 

20 

-0.0472 

0.0138 

0.0750 

0.1368 

30 

-0.0965 

-0.0353 

0.0267 

0.0885 

Bottom  End 

Positive  Diagonal 

0 

10 

20 

30 

Negative 

Diagonal 

0 

-0.0662 

-0.0713 

-0.1053 

-0.1428 

10 

-0.0072 

-0.0383 

-0.0756 

-0.1133 

20 

0.0289 

-0.0085 

-0.0460 

-0.0838 

30 

0.0592 

0.0216 

-0.0164 

-0.0543 
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Table  15.  Lateral  displacement  results  for  the  closing  cases. 


Positive  Diagonal 

0 

10 

20 

30 

Top  End 

Negative 

Diagonal 

0 

-0.2269 

-0.1217 

-0.0155 

0.0596 

10 

-0.2258 

-0.1170 

-0.0519 

0.0140 

20 

-0.2236 

-0.1614 

-0.0997 

-0.0384 

30 

-0.2712 

-0.2096 

-0.1480 

-0.0863 

Positive  Diagonal 

0 

10 

20 

30 

Bottom  End 

Negative 

Diagonal 

0 

0.1391 

0.0746 

0.0095 

-0.0365 

10 

0.1384 

0.0717 

0.0318 

-0.0086 

20 

0.1371 

0.0989 

0.0611 

0.0235 

30 

0.1662 

0.1285 

0.0907 

0.0529 

Figure  40  presents  a  typical  curve  of  the  miter  end  behavior  for  a  positive 
diagonal  prestressing  of  20  ksi  and  various  negative  diagonal  prestressing 
when  the  gate  is  opening  and  closing.  The  plot  demonstrates  the  relation 
between  the  diagonal  prestressing  and  the  opening  and  closing  behavior  of 
the  gate.  As  mentioned  in  the  previous  parametric  cases,  the  positive  diag¬ 
onal  prestressing  forces  the  top  portion  of  the  miter  end  to  move  in  the 
upstream  direction,  and  vise  versa  when  the  prestressing  increases  in  the 
negative  diagonal. 

That  implies  that  an  increment  in  the  positive  diagonal  prestressing  pro¬ 
duces  an  upstream  movement  of  the  top  portion  of  the  gate.  This  move¬ 
ment  then  produces  a  better  behavior  of  the  gate  when  closing  because 
when  the  gate  is  closing,  the  top  of  the  miter  gate  moves  downstream.  On 
the  other  hand,  an  increment  in  the  negative  diagonal  prestressing  results 
in  a  better  behavior  of  the  gate  when  opening  because  the  lateral  displace¬ 
ment  produced  by  the  prestressing  in  the  negative  diagonal  is  reduced 
when  the  gate  opens.  This  means  that  the  gate  will  experience  a  relative 
lateral  displacement  when  opening  and  closing,  but  an  adequate  diagonal 
prestressing  combination  can  maintain  the  lateral  displacements  between 
tolerable  limits. 
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Open  and  Close  Lateral  Displacement 
Pos.  Diag.  20  ksi  and  Neg.  Diag.  Vary 


Downstream  Lateral  Displacement  [ft]  Upstream 


Open  20  0  Close  20  0  Open20_10  Close  20_  10 

-A-  Open  20  20  Close  20_20  Open  20_30  Close  20_30 

Figure  40.  Typical  miter  end  behavior. 

Figure  41  summarizes  all  the  lateral  displacements  for  the  positive  and 
negative  diagonal  combinations.  In  order  to  obtain  the  optimal  diagonal 
prestressing  combination  when  the  gate  opens  and  closes,  a  maximum 
lateral  displacement  limit  had  to  be  defined  (tolerance).  For  example,  if 
this  maximum  limit  is  1.0  in.  of  lateral  displacement  in  the  miter  end,  then 
based  on  Figure  41  any  diagonal  prestressing  combination  will  satisfy  the 
maximum  displacements  requirement  for  opening  and  closing  the  gate. 

Figures  42  and  43  correspond  to  the  opening  and  closing  cases,  respect¬ 
ively.  The  first  step  is  to  set  the  maximum  lateral  displacement  limits  in 
the  plot;  second;  select  the  diagonal  prestressing  combinations  that  satisfy 
the  maximum  lateral  displacements  limits;  and  third,  choose  all  the 
prestressing  combinations  that  satisfy  the  predefined  tolerance. 

Figure  44  presents  the  relative  lateral  displacement  of  the  miter  end  when 
the  negative  diagonal  prestressing  varies  and  the  positive  diagonal 
prestressing  is  maintained  constant  in  each  series. 
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Movement  +  Diagonal  Prestressing 
Lateral  Displacement  [U2G] 
[Open  and  Close] 


♦  0  ksi  (Top)_0 

-■-10  ksi  (Top)  O 

20  ksi  (Top)_0 

30  ksi  (Top)  O 

•  0  ksi  (Bot)  O 

—^10  ksi  (Bot)  O 

20  ksi  (Bot)  O 

30  ksi  (Bot)  O 

-*-0ksi(Top)_C 

10  ksi  (Top)  C 

-■—20  ksi  (Top)_C 

A  30  ksi  (Top)_C 

X  0  ksi  (Bot)C 

10  ksi  (Bot)  C 

•  20  ksi  (Bot)  C 

■  30  ksi  (Bot)C 

Figure  41.  Lateral  displacements  with  different  diagonal  prestressing 

combination. 


Movement  +  Diagonal  Prestressing 
Lateral  Displacement  [U2G] 
[Open] 


Displacement  [ft] 

♦  0  ksi  (Top)  ■  10  ksi  (Top)  20  ksi  (Top)  30  ksi  (Top) 

T  0  ksi  (Bot)  1  10  ksi  (Bot)  20  ksi  (Bot)  30  ksi  (Bot) 

Figure  42.  Lateral  displacement-  different  diagonal  prestressing  combination, 

gate  open. 
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Movement  +  Diagonal  Prestressing 
Lateral  Displacement  [U2G] 
[Close] 


Displacement  [ft] 

0  ksi  (Top)_C  1 0  ksi  (Top)_C  20  ksi  (Top)_C  -±-  30  ksi  (Top)_C 

-*-0ksi(Bot)_C  -*-10ksi(Bot)_C  -#-  20  ksi  (Bot)_C  — 1—30  ksi(Bot)_C 


Figure  43.  Lateral  displacement-different  diagonal  prestressing  combination, 

gate  closed. 


Operation  +  Water  Resistance 
Relative  Lateral  Displacement 


Relative  Lateral  Displacement  [ft] 


•  10  ksi  Open 

•  30  ksi  Open 


■  10  ksi  Close 

■  30  ksi  Close 


♦  0  ksi  Open 
)K  20  ksi  Open 


"  0  ksi  Close 
20  ksi  Close 


Figure  44.  Relative  lateral  displacement  (negative  diagonals  vary). 
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Figure  45  presents  the  relative  lateral  displacement  when  the  positive 
diagonal  prestressing  varies  and  the  negative  diagonal  prestressing  is 
maintained  constant  in  each  series.  These  plots  can  be  used  to  determine 
the  zero  relative  displacement  for  a  given  diagonal  prestressing  when  the 
gate  is  opening  and  closing. 


Operation  +  Water  Resistance 
Relative  Lateral  Displacement 


♦  0  ksi  Open  ■  0  ksi  Close  *  10  ksi  Open  X  10  ksi  Close 

)K  20  ksi  Open  •  20  ksi  Close  I  30  ksi  Open  ■  30  ksi  Close 

Figure  45.  Relative  lateral  displacement  (positive  diagonals  vary). 


5.5.5  Case  5:  Quoin  block  analysis 

Quoin  block  will  deteriorate  during  the  service  live  of  the  structure.  The 
development  of  the  deterioration  will  cause  that  the  axial  loads  transfer  to 
the  wall  be  diminished  causing  an  increase  of  vertical  reactions  carry  by 
the  pintle  region.  For  the  analysis  a  constant  value  of  diagonal  prestressing 
was  used.  The  pre-stressing  in  the  diagonals  were  20  ksi  for  the  positive 
diagonal  and  25  ksi  for  the  negative  diagonal. 

Figure  46  shows  a  schematic  figure  for  this  case,  and  the  assumed 
boundary  conditions  for  this  case  are  shown  in  Table  16.  The  different 
cases  used  to  simulate  the  quoin  block  deterioration  are  presented  and 
described  in  Table  17. 
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Gate 

Component 

U1 

U2 

U3 

UR1 

UR2 

UR3 

CSYS 

Gudgeon  pin 

0 

0 

— 

— 

— 

— 

Global 

Pintle 

0 

0 

0 

0 

0 

— 

Global 

Quoin 

0 

0 

— - 

— 

— 

— 

Local 

Miter  end 

0 

0 

— - 

— 

— 

— 

Local 
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Table  17.  Quoin  block  boundary  condition  variation. 


BC 


Analysis  Description 


Miter 

Gate 


Q100: 

Full  quoin  boundary  condition. 


Miter 

Gate 


Q75T: 

75  %  of  the  quoin  area  measured  from  top  to  bottom. 


Miter 

Gate 


Q50T: 

50  %  of  the  quoin  area  measured  from  top  to  bottom. 


Miter 

Gate 


Q25T: 

25  %  of  the  quoin  area  measured  from  top  to  bottom. 


Miter 

Gate 


Q75B: 

75  %  of  the  quoin  area  measured  from  bottom  to  top. 


Miter 

Gate 


Q50B: 

50  %  of  the  quoin  area  measured  from  bottom  to  top. 


Miter 

Gate 


Q25B: 

25  %  of  the  quoin  area  measured  from  bottom  to  top. 


Miter 

Gate 


Q50C: 

50  %  of  the  quoin  area  measured  from  the  center. 


Miter 

Gate 


Q25C: 

25  %  of  the  quoin  area  measured  from  the  center. 


Miter 

Gate 


QO: 

No  quoin  boundary  condition.  Free  Nodes. 


When  the  gate  is  closed,  the  water  level  in  the  upstream  face  of  the  gate  is 
the  maximum  head  and  the  water  level  in  the  downstream  face  of  the  gate 
increases  to  reach  the  upstream  level.  A  total  of  ten  steps  were  used  to 
represent  the  filling  of  the  lock  chamber.  Each  step  represents  an  increase 
in  the  water  level  of  the  downstream  side  of  the  gate.  In  the  final  step  the 
hydrostatic  pressure  in  the  upstream  surface  of  the  gate  is  balanced  with 
the  downstream  water  pressure.  Figure  47  presents  the  different  water 
levels  used  in  the  analysis.  After  these  ten  steps  the  stresses  in  the  pintle 
zone  and  in  the  critical  members  were  investigated. 


ERDC/ITL  TR-09-1 


62 


5.6  Stress  concentration  at  pintle  zone 

Table  18  summarizes  the  stresses  in  the  pintle  zone  for  each  of  the  cases. 
The  results  reveal  an  increase  in  the  stresses  in  the  pintle  zone  for  the 
cases  where  the  quoin  block  is  deteriorated  from  top  to  bottom  (Q#T)  and 
from  the  center  (Q#C)  compared  with  Qioo,  which  represents  a  full  action 
of  the  quoin  block.  Figure  48  presents  the  results  obtained  for  quoin 
blocks  deteriorating  from  top  to  bottom.  This  plot  demonstrates  clearly  an 
increment  in  the  stresses  at  the  pintle  zone  when  the  quoin  block  area  that 
transfers  the  axial  load  is  reduced.  Figure  49  presents  the  results  for  the 
quoin  block  cases  deteriorating  from  the  center  up  and  down.  When  the 
quoin  block  starts  to  be  ineffective,  the  stresses  in  the  pintle  zone  begin  to 
rise  as  a  function  of  the  effective  area.  The  maximum  stress  value  was 
obtained  in  the  analysis  when  the  quoin  block  was  inactive  in  the  structure 
(Qo).  Some  of  the  stresses  obtained  in  the  different  analyses  were  above 
the  yield  strength  of  the  members  (50  ksi).  The  cases  when  the  quoin  block 
is  deteriorating  from  the  bottom  portion  to  the  top  (Q#B)  do  not  produce 
an  increment  in  the  stresses  at  the  pintle  zone.  The  values  obtained  in  the 
results  are  around  10  ksi  (Figure  50).  After  a  careful  study  of  the  structure 
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behavior  in  the  model,  a  concentration  of  stresses  in  the  quoin  block 
effective  area  was  detected.  This  effect  produces  increments  in  the  stresses 
of  the  structural  elements  in  the  gate.  Table  19  and  Figure  51  present  the 
stresses  in  the  critical  zone  for  these  cases.  The  stress  values  are  closer  to 
the  yield  strength  of  the  members,  and  a  quoin  block  effective  area  less 
than  25  percent  can  exceed  the  member  capacity.  Figure  52  demonstrates 
the  stress  concentration  in  the  quoin  load  transfer  area  for  Case  Q25B. 


Table  18.  Summary  stresses  at  pintle  zone. 


Quoin 

Case 

Stress 

ksi 

Quoin 

Case 

Stress 

ksi 

Q100 

10.7 

Q25C 

38.1 

Q75T 

29.3 

Q75B 

10.8 

Q50T 

47.0 

Q50B 

9.4 

Q25T 

63.8 

Q25B 

9.8 

Q50C 

29.5 

Q0 

80.2 

Stress  at  Pintle  Zone 


Quoin  Case 


Figure  48.  Stress  concentration  in  the  pintle  zone  for  cases  from  top  to  bottom. 
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Stress  at  Pintle  Zone 
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Q50C  Q25C 
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Figure  49.  Stress  concentration  in  the  pintle  zone  for  cases  from  center 

up  and  down. 
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Figure  50.  Stress  concentration  in  the  pintle  zone  for  cases  from  bottom  to  top. 
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Table  19.  Stresses  at  the  critical  zone. 


Quoin  Case 

Stress  ksi 

Q100 

10.7 

Q75B 

10.8 

Q50B 

18.4 

Q25B 

45.6 

QO 

80.2 

Stress  at  Critical  Zone 


Quoin  Case 


Figure  51.  Stress  concentration  in  the  critical  zone  for  cases  from  bottom  to  top. 


The  results  obtained  demonstrate  that  deterioration  in  the  quoin  block 
impeding  the  transfers  of  axial  loads  to  the  wall  can  affect  drastically  the 
state  of  stresses  in  the  pintle  region.  If  the  deterioration  is  severe,  the 
stresses  can  reach  undesirable  stress  levels.  The  location  of  maximum 
stress  concentration  depends  on  the  quoin  deterioration  area:  if  the 
deterioration  occurs  in  the  pintle  area  (top  -  bottom),  the  maximum 
stresses  will  be  produced  in  the  pintle  zone;  and  if  the  deterioration  is  in 
the  gudgeon  pin  (bottom  -  top),  the  maximum  stress  will  be  produced  in 
the  elements  near  the  quoin  block  effective  area  end. 
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Figure  52.  Stress  concentration  in  the  critical  zone  for  Case  Q25B. 
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6  Conclusions  and  Recommendations 

6.1  3-D  finite  element  analysis  of  miter  gate  with  multiple  cracks 

This  Technical  Report  has  described  the  capability  of  calculating  stress 
intensity  factors  using  state-of-the-art  analytical  techniques.  The  calcu¬ 
lation  of  the  stress  intensity  factors  can  be  obtained  from  a  3-D  finite 
element  analysis.  3-D  finite  element  meshes  can  be  automatically 
generated  within  commercially  available  applications  using  IGES- 
formatted  CAD  drawings. 

The  first  example  showed  that  the  J-integral  approach  successfully  pre¬ 
dicted  the  stress  intensity  factors  for  a  semi-infinite  plate  with  through¬ 
thickness  crack,  double-edge  cracks,  single-edge  crack,  and  cracks  growing 
from  round  holes.  The  example  showed  a  difference  of  less  than  2  percent 
between  the  numerically  calculated  values  and  the  ones  calculated  with 
closed-form  solutions. 

In  the  second  example,  the  3-D  meshing  technique  utilizing  CAD  drawings 
imported  into  the  finite  element  application  was  demonstrated.  The 
meshing  procedure  allows  the  gates  with  all  its  components  to  be  modeled 
completely  without  needing  to  make  any  simplifying  assumptions  regard¬ 
ing  the  geometry.  Once  again,  the  stress  intensity  factors  of  the  upstream 
flanges  of  a  horizontal  girder  and  a  vertical  diaphragm  were  calculated. 

For  this  example,  the  stress  intensity  factors  showed  that  the  cracks  will 
not  grow.  Therefore  repairs  can  be  postponed. 

The  J-integral  methodology  can  certainly  be  used  to  evaluate  cracked 
hydraulic  steel  structures  and  to  determine  remaining  life  based  on  the 
operation  conditions.  The  methodology  can  also  be  used  to  decide  the 
criticality  of  the  cracking  problem  and  whether  repairs  are  needed. 
Suggested  repairs  can  also  be  evaluated  numerically  to  predict  reliability 
and  safety  conditions.  These  techniques  when  coupled  with  a  library  of  all 
of  the  Corps  hydraulic  structures  will  allow  quick  analysis  whenever 
problems  are  found. 

Further  studies  regarding  crack  initiation  should  be  performed.  Crack 
initiation  in  welded  steel  structures  is  influenced  by  fatigue  loading, 
connection  details,  thermal  residual  stresses,  and  material  properties. 
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Fracture  modes  II  and  III  should  also  be  evaluated  to  establish  the  contri¬ 
bution  of  shear  and  out-of-plane  shear  to  the  stress  intensity  factor. 

6.2  Numerical  investigation  of  diagonals  prestress 

The  comprehensive  miter  gate  diagonals  prestress  analysis  demonstrated 
the  importance  of  maintaining  the  adequate  levels  of  prestress.  Opening 
and  closing  the  gate  produces  an  additional  lateral  displacement.  Relax¬ 
ation  will  develop  over  time  and  vertical  and  lateral  displacement  on  the 
miter  end  will  increase,  causing  misalignment  and  additional  torsional 
stresses  on  the  pintle  region. 

Optimum  values  of  prestressing  can  easily  be  calculated  using  numerical 
simulations.  Once  the  different  diagonal  prestressing  models  are  obtained, 
a  series  of  curves  relating  prestressing  levels  with  relative  displacements 
can  be  accomplished,  then  for  a  predefined  tolerance  the  optimum 
prestressing  levels  can  be  easily  extracted  from  these  graphs.  During  the 
gravity  analysis  the  results  revealed  a  twist  in  the  miter  gate  end.  The 
quoin  end  of  a  miter  gate  is  held  vertical  by  the  pintle  and  gudgeon  pin, 
leaving  the  miter  end  free  to  twist  out  of  vertical  alignment.  This  twist  was 
produced  because  the  center  of  gravity  and  the  shear  center  are  not  in  the 
same  vertical  plane. 

The  case  of  gravity  with  positive  diagonal  showed  that  the  greater  the 
prestressing  in  the  positive  diagonals,  the  smallest  the  vertical  deflection 
becomes.  Zero  or  positive  vertical  deflections  can  be  accomplished  for 
prestress  values  higher  than  the  yield  strength;  therefore,  a  negative  diag¬ 
onal  is  required  to  provide  additional  vertical  stiffness.  Prestressing  of  the 
positive  diagonal  contributes  to  reduce  the  twisting  generated  by  the 
gravity  load. 

The  gravity  and  positive  and  negative  diagonal  case  shows  that  the  greater 
the  prestressing  in  the  negative  diagonals,  the  greater  the  vertical  displace¬ 
ment  becomes.  With  this  analysis  the  optimum  prestressing  levels  can  be 
obtained  by  using  plots  generated  with  different  levels  of  prestressing. 

The  case  involving  gravity  with  positive  and  negative  diagonals  during  gate 
operation  showed  that  the  gate  will  experience  a  relative  lateral  displace¬ 
ment  when  opening  and  closing,  but  an  adequate  diagonal  prestressing 
combination  can  maintain  the  lateral  displacements  between  tolerable 
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limits.  These  limits  can  be  found  from  the  figures  included  showing  the 
vertical  and  lateral  displacements  fro  different  diagonal  prestressing. 

Quoin  block  deterioration  analysis  demonstrated  that  deterioration  in  the 
quoin  block  can  affect  drastically  the  state  of  stresses  in  the  pintle  region. 
If  the  deterioration  is  severe,  the  stresses  can  reach  undesirable  stress 
levels.  The  location  of  the  stress  concentrations  depends  on  the  quoin 
deterioration  area:  if  the  deterioration  occurs  in  the  pintle  area  (top  - 
bottom),  the  maximum  stresses  will  be  generated  in  the  pintle  zone;  and  if 
the  deterioration  is  in  the  gudgeon  pin  (bottom  -  top),  the  maximum 
stress  will  be  produced  in  the  elements  near  the  quoin  block  effective  area 
end. 

6.3  Recommendations 

A  comprehensive  3-D  finite  element  analysis  can  be  used  to  determine  the 
stress  intensity  factor  for  cracked  hydraulic  steel  structures.  This  factor  is 
essential  for  the  assessment  and  determination  of  the  remaining  life  of  a 
fatigue-induced  cracking  problem. 

The  advanced  numerical  modeling  can  also  be  used  to  asses  the  optimum 
values  of  diagonal  prestressing  for  new  and  existing  miter  gates,  elim¬ 
inating  the  possibility  of  diagonal  connection  failures. 

It  is  essential  that  most  of  the  Corps  of  Engineers  navigation  hydraulic 
steel  structures  be  evaluated  numerically.  This  type  of  analysis  will  help  in 
the  maintenance  and  evaluation  of  repairs  when  necessary.  Furthermore, 
considering  the  aging  navigation  infrastructure  in  this  country,  by  con¬ 
ducting  these  advance  numerical  evaluations,  the  Corps  will  be  proactive 
instead  of  reactive  to  possible  problems.  By  being  proactive,  the  Corps  will 
obtain  substantial  savings  when  performing  repairs  and  maintenance,  or 
when  a  critical  problem  arises. 
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